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16.  Ak»tr««r 

^  model  developed  previously  as  part  of  the  Hazard  Assessment  Computer  System  to 
predict  the  discharge  race  of  a  chemical  cargo  from  a  punctured  ship  tank  has 
been  revised  and  verified.  The  computerized  model  can  now  treat:  (1)  a  puncture 
above,  below,  or  intersecting  the  waterline;  (2)  vacuum  relief  valve  actuations; 
(3)  simultaneous  liquid  and  gas  discharge;  (4)  volatile  or  nonvolatile  cargo;  and 
(5)  air  or  water  ingestion  into  Che  tank,  to  permit  additional  cargo  to  be  dis¬ 
charged.  The  model  is  based  on  an  analysis  of  such  thermodynamic  and  fluid  flow 
phenomena  as  two-phase  and  choked  flow  and  realistic  evaporation  rates  within  the 
tank.  The  Chemical  Properties  File  was  examined  for  its  applicability  to  the 
revised  model,  and  some  desirable  improvements  were  indicated.  Experimental  veri¬ 
fication  of  the  model  consisted  of  a  series  of  extensively  instrumented  tests  of 
small  tanks  to  investigate:  (1)  discharge  coefficients  of  irregularly  shaped 
punctures;  (2)  air  and  water  ingestion;  and  (3)  discharge  of  volatile  and  boiling 
liquids.  The  good  comparison  of  Che  model  to  the  tests  adequately  verified  the 
model . i . 
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SUMMARY 


? 


This  final  report  covers  all  four  tasks  of  a  project  to 
revise  and  verify  experimentally  the  Venting  Rate  Model  of  the 
Hazards  Assessment  Computer  System.  The  report  documents  (1) 
all  experimental  procedures  and  representative  test  data,  and 
(2)  th-\  analysis,  development,  and  verification  of  the  final 
form  of  the  model. 

Background .  An  analytical/computer  model  had  been  previ¬ 
ously  developed  for  the  U.  S.  Coast  Guard  for  use  in  predict¬ 
ing  the  discharge  of  a  liquid  or  gas  cargo  from  a  punctured 
ship  tank.  A  later  independent  review  found  that  the  model 
contained  a  number  of  errors  and  inappropriate  assumptions. 

The  review  generally  concluded  that  the  model  should  be  re¬ 
vised  and  verified  experimentally.  The  present  program  was 
therefore  sponsored  by  the  Coast  Guard  to  accomplish  the  rec¬ 
ommended  revision  and  verification. 

Tasks  1  and  2  -  Review  and  Reformulation  of  the  Model. 
During  the  first  two  tasks  of  this  program,  the  Venting  Rate 
Model  and  other  relevant  literature  were  reviewed,  and  the 
model  was  then  completely  reformulated.  The  new  model  is  ap¬ 
plicable  to  nearly  every  practical  combination  of  puncture 
location  ,  cargo  properties,  intermittent  outflows,  and  relief 
valve  actuations.  In  summary,  the  new  model  incorporates: 

o  puncture  locations  above  or  below  the  waterline? 
o  vacuum  relief  valve  actuations; 

o  possibility  of  air  in  the  cargo  vapor  space; 
o  changes  of  phase  during  the  discharge  of  a 
volatile  cargo? 

o  choked  flow  discharges  of  gases  and  volatile  liquids; 
o  simultaneous  liquid  and  gas  discharges; 

o  air  or  water  ingestion  through  the  puncture;  and 
o  effects  of  realistic  evaporation  rates  of  volatile 
liquids  within  the  tank. 


Temperature  stratification  of  a  volatile  liquid,  one  of  the 
phenomena  suggested  by  the  review  as  perhaps  being  important 
enough  to  be  included  in  the  model,  was  investigated  experi¬ 
mentally  and  found  to  be  negligible.  Thus,  stratification  is 
not  included  in  the  new  model. 

The  model  has  been  programmed  for  computerized  solution. 
Program  listings  and  flow  charts  are  included  in  this  report. 

The  Chemical  Properties  File  was  also  reviewed  for  its 
applicability  to  the  revised  Venting  Rate  Model.  Data  corre¬ 
lations  for  cargo  properties  of  nonvolatile  cargos  were  found 
to  be  satisfactory.  However,  several  serious  limitations  were 
discovered  for  volatile  cargos;  namely;  (1)  the  compressibility 
and  the  specific  heat  of  the  saturated  vapor  should  be  in¬ 
cluded;  (2)  the  vapor  pressure  vs  temperature  correlations  are 
moderately  inaccurate;  and  (3)  a  correlation  giving  the  latent 
heat  of  vaporation  as  a  function  of  temperature  should  be  in¬ 
cluded.  More  accurate  and  complete  correlations  for  six  rep¬ 
resentative  cargos  were  also  developed. 

Task  3  -  Experimental  Program.  The  third  task  of  the 
program  was  an  extensive  series  of  model-scale  experiments. 

The  overall  objectives  of  these  tests  were: 

o  obtain  data  to  evaluate  the  assumptions  used  in 
developing  the  analytical  models; 
o  gain  insight  into  air  and  water  ingestion;  and 

o  acquire  discharge  rate  data  for  a  variety  of  test 

conditions  to  validate  the  overall  Venting  Rate 
Model. 

The  first  test  series  determined  the  discharge  coeffi¬ 
cients  for  a  variety  of  realistically  shaped  punctures  over  a 
range  of  Reynolds  numbers.  These  tests  showed  that  a  discharge 
coefficient  of  0.65  is  adequate  for  nearly  all  puncture  shapes; 
the  exceptions  are  vertical  or  horizontal  slot-like  punctures 
with  ragged  edges  pointing  outwards,  for  which  the  recommended 
discharge  coefficient  is  0.825. 


The  second  test  series  was  designed  to  study  the  inges¬ 
tion  of  air  or  water  (for  submerged  punctures)  through  the 
puncture.  It  was  found  that  air  or  water  ingestion  can  result 
in  the  discharge  of  more  cargo  than  would  be  predicted  by  the 
previous  Venting  Rate  Model  when  the  ambient  pressure  exceeds 
the  tank  pressure  at  the  puncture  location.  The  tests  were 
used  to  obtain  physical  insight  about  air  and  water  ingestion, 
from  which  relatively  simple  but  realistic  analytical  models 
were  developed. 

The  third  test  series  studied  the  discharge  of  volatile 
cargos.  The  specific  objectives  of  these  tests  were:  (1) 
evaluate  the  importance  of  liquid  stratification;  (2)  gain  in¬ 
sight  and  obtain  quantitative  data  about  evaporation  within 
the  tank;  (3)  evaluate  the  effects  of  phase  changes  of  the 
discharge  (two-phase  outflow) ;  and  (4)  obtain  discharge  data 
(e.g.,  time  durations,  temperatures,  tank  pressure,  mass  rate 
of  outflow)  which  could  be  used  for  the  overall  verification 
of  the  revised  analytical/computer  model. 

As  indicated  above,  these  tests  showed  that  liquid  strati¬ 
fication  was  negligible.  The  effects  of  evaporation  rate,  how¬ 
ever,  were  found  to  be  significant  when  the  vacuum  relief  valve 
is  jammed  shut  or  when  the  cargo  is  transported  at  super-atmo¬ 
spheric  pressure.  Evaporation  within  the  tank  tends  to  keep 
the  partial  pressure  of  the  vapor  in  the  tank  at  some  value 
intermediate  between  the  saturation  pressure  of  the  liquid  and 
the  pressure  that  corresponds  to  the  expansion  of  a  gas  whose 
mass  does  not  increase  as  the  vapor  space  volume  increases 
during  the  discharge.  Since  the  tank  pressure  is  an  important 
factor  in  predicting  the  discharge  rate,  a  new  model  of  evapo¬ 
ration  rates  appropriate  for  closed  tanks  was  developed.  The 
tests  also  showed  that  the  discharge  rate  of  liquids  was  lim¬ 
ited  by  changes  of  phase  to  smaller  values  than  would  be  pre¬ 
dicted  by  the  single-phase  flow  relations  used  in  the  previous 
Venting  Rate  Model.  Two-phase  flow  models  were  therefore  de¬ 
veloped. 


Task  4  -  Model  Verification.  In  this  task,  the  Venting 
Rate  Model  was  compared  to  the  test  data  and  revised  as  indi¬ 
cated  by  the  comparisons.  This  was,  in  fact,  a  continuing 
process  during  the  entire  project. 

Each  submodel  of  the  overall  model  was  first  compared  to 
representative  test  data  relevant  to  that  model  and  then  cor¬ 
rected  when  necessary.  In  addition,  the  tests  were  used  to 
determine  the  best  values  for  any  empirical  constants  appear¬ 
ing  in  the  model.  It  should  be  noted  that  no  submodel  re¬ 
quired  more  than  one  such  empirical  constant  and  that  the 
values  of  the  constants  were  all  physically  reasonable. 

The  entire  model  was  then  compared  to  tests  for  many  dif¬ 
ferent  kinds  of  cargo  discharge: 

o  nonvolatile  liquids,  with  and  without  operable 
vacuum  relief  valves; 

o  air  ingestion  during  the  discharge  of  volatile 
and  nonvolatile  liquids; 

o  water  ingestion  during  the  discharge  of  nonvolatile 
liquids  heavier  and  lighter  than  water;  and 

o  volatile  liquids  transported  at  slightly  super- 
atmospheric  pressures. 

The  comparisons  of  the  model  to  these  tests  were  very  good. 

It  is  concluded  that  the  revised  model  is  now  satisfactory 
for  use  in  the  Hazards  Assessment  Computer  System. 
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I .  INTRODUCTION 


As  part  of  the  Hazards  Assessment  Computer  System  (HACS) 
of  the  Chemical  Hazards  Response  Information  System,  a  model 
was  previously  developed  to  predict  the  rate  of  liquid  or  gas¬ 
eous  cargo  discharged  to  the  environment  by  an  accidental  rup¬ 
ture  of  a  ship  tank  [1] .  An  independent  review  of  this  model 
later  revealed  several  errors  in  its  basic  equations  and  a 
lack  of  realism  in  some  of  its  underlying  assumptions  [2] . 

The  Coast  Guard  has  sponsored  the  present  program  to  correct 
these  deficiencies  and  to  validate  the  corrected  model  experi¬ 
mentally. 

Errors  in  the  basic  equations  and  incorrect  assumptions 
made  in  developing  the  previous  model  included: 

o  the  thermodynamic  energy  balance  used  to  predict 
tank  pressure  and  temperature  for  volatile  cargos 
was  incorrect; 

o  incompressible  flow  equations  were  used  to  predict 
the  venting  of  a  volatile  liquid,  even  though  the 
liquid  partially  vaporizes  during  the  venting; 
o  the  discharge  was  assumed  to  cease  when  the  pres¬ 
sure  differential  across  the  puncture  was  zero, 
even  if  the  liquid  cargo  level  were  still  above 
the  puncture;  that  is,  ingestion  of  air  or  water 
through  the  puncture  was  ignored  although  this 
would  allow  the  remaining  liquid  cargo  above  the 
puncture  to  be  discharged;  and 
o  incorrect  values  were  recommended  for  the  discharge 
coefficients  of  the  puncture. 

Other  stated  or  implied  assumptions  that  may  be  unrealistic  or 
inappropriate  included: 

o  no  provision  was  made  for  vacuum  relief  valves; 
o  all  punctures  were  above  the  waterline; 


o 


simultaneous  venting  of  liquid  and  gas  did  not 
occur; 

o  air  was  not  present  in  the  vapor  space; 

o  the  specific  heat  of  superheated  vapor  was  used 

to  compute  enthalpy  changes  of  saturated  vapor; 

o  evaporation  at  the  liquid-gas  interface  in  the 

tank  was  computed  by  assuming  the  vapor  remained 
saturated;  and 

o  temperature  stratification  in  the  liquid  resulting 
from  evaporation  was  negligible. 

The  consequences  of  these  errors  and  assumptions  could  not 
always  be  determined  solely  by  examining  the  model  or  its  pre 
dictions.  For  example,  potential  errors  in  the  evaporation 
model  could  not  be  evaluated  because  a  more  realistic  model 
was  not  available.  Further,  the  overall  accuracy  of  the  vent 
ing  rate  model  was  unknown. 

The  present  program  was  designed  to  correct  these  errors 
and  to  validate  the  models.  It  had  the  following  scope. 

1.  Review  the  venting  rate  model  and  any  other  rele¬ 
vant  literature,  with  the  aim  of  uncovering  errors 
inconsistencies,  and  unrealistic  assumptions. 

2.  Correct  and  reformulate  the  venting  rate  model  and 
reprogram  it  for  computerized  solution. 

3.  Review  the  Chemical  Properties  File  of  HACS  for 
applicability  to  the  corrected  venting  rate  model. 

4.  Conduct  a  sensitivity  analysis  of  the  computed  re¬ 
sults  to  the  input  variables. 

5.  Design  and  conduct  a  model-scale  test  program  to 
obtain  venting  rate  data  for  tank  conditions  and 
cargo  conditions  analogous  to  full  scale. 

Revise  the  venting  rate  model  as  indicated  by  the 
correlation  of  the  test  data  and  the  model  predic¬ 
tions. 
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These  efforts  were  arranged  into  four  tasks.  Task  1  covered 
the  first  effort  given  above?  Task  2  covered  the  next  three; 
Task  3  covered  the  fifth  effort;  and  Task  4  covered  the  last 
effort.  To  present  a  more  readable  narrative,  this  final  re¬ 
port  is  not  arranged  strictly  by  tasks.  The  complete  analyti¬ 
cal  model  development  is  presented  first  (Tasks  1,  2,  and  4) . 
Then  the  test  program  (Task  3)  is  described.  Comparisons  of 
the  tests  and  the  model  predictions  are  given  next,  followed 
by  the  conclusions  and  appendices,  which  incorporate  the  sen¬ 
sitivity  analysis  and  the  review  of  the  Chemical  Properties 
File . 
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II.  TASK  II.  REFORMULATION  OF  MODEL 


The  venting  rate  model  was  reformulated  in  detail  in  the 
Interim  Report  [3] .  A  summary  of  that  development  is  presented 
in  this  section,  with  a  complete  derivation  presented  in  Appen¬ 
dix  A.  More  complete  details  of  the  air  and  water  ingestion 
submodels  and  the  evaporation  rate  model  will  be  given,  however, 
since  they  were  not  fully  developed  at  the  time  the  Interim 
Report  was  issued. 

II.l  Reformulation  of  Tank  Energy  Balance 

The  total  mass  of  cargo  discharged  can  sometimes  be  easily 
estimated  (for  example,  the  quantity  of  liquid  above  the  punc¬ 
ture)  .  But  the  venting  rates  and  the  time  required  for  the 
venting  can  only  be  determined  by  coupling  energy  and  mass 
balances  which  describe  the  time-varying  pressure  and  tempera¬ 
ture  of  the  cargo  within  the  tank  cargo  to  models  of  the  dis¬ 
charge  rate  through  the  puncture.  As  in  AMSHAH  [1],  two  limit¬ 
ing  conditions  are  assumed  in  formulating  the  energy  balances: 

(1)  all  processes  occurring  in  the  tank  are  adiabatic, 
or 

(2)  all  processes  occurring  in  the  tank  are  isothermal. 

These  assumptions  eliminate  the  need  to  consider  energy  trans¬ 
fers  with  the  outside  and  thus  greatly  simplify  the  models. 
Furthermore,  predictions  made  using  these  assumptions  will 
bound  the  actual  tank  conditions.  In  addition,  isothermal  and 
adiabatic  conditions  are  practically  identical  for  a  nonvola¬ 
tile  cargo,  since  the  temperature  changes  little  during  the 
venting,  and  so  for  this  case  the  assumption  of  which  process 
actually  occurs  is  not  important. 

The  critical  review  [2]  of  AMSHAH 's  model  concluded  that 
temperature  stratification  as  a  result  of  heat  and  mass  trans¬ 
fer  between  the  liquid  and  vapor  phases  of  the  cargo  within 
the  tank  might  need  to  be  accounted  for  in  the  energy  balances. 
Experimental  results  to  be  described  later  show,  however,  that 
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the  temperature  stratification  is  small  and  has  only  a  negli¬ 
gible  effect  on  the  tank  energy  balance,  even  for  very  rapid 
venting  of  volatile  cargos.  In  agreement  with  the  original 
models,  then,  a  "well  stirred  reactor"  approach  is  assumed. 

Considering  the  control  volume  shown  in  Figure  II. 1  and 
assuming  for  the  moment  that  the  processes  within  the  tank 
are  adiabatic,  the  correct  form  of  the  energy  balance  for  the 
cargo  in  the  tank  is  [31 : 


<Vpl  +  «V=PV  *  maSa>  33 
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(II. 1) 


The  negligible  changes  in  gravitational  potential  energy  and 
kinetic  energy  have  been  ignored  in  this  equation  [1,2,3]. 

The  bar  over  Cpv  is  meant  to  emphasize  that  the  enthalpy 
change  of  saturated  vapor  (dhy  *  CpvdT)  is  the  relevant  quan¬ 
tity  for  the  vapor,  for  which  Cp  of  superheated  vapor  is  a 
poor  approximation.  For  isothermal  processes,  equation  (II.l) 
reduces  simply  to  T  =  constant  »  initial  temperature. 

The  mass  balances  which  must  be  used  with  equation  (II.l) 
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Since  VT  *  v^M^  +  vviV'  where  the  specific  volume  of  the  va¬ 
por  is  based  on  its  partial  pressure  in  the  air/vapor  mixture 
in  the  vapor  space,  and  since  dv^/dt  *  0,  equations  (II. 2a) 
and  (II. 2b)  can  be  combined  to  give  a  relation  for  the  evapo¬ 
ration  rate: 


NOMENCLATURE  FOR  TANK  ENERGY  BALANCE 
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r  dv  dv  — J  / 

We  -  [V't.o  *  Wo  -  mL  TIT  -  «v  *3rj/(vv  *  VLJ  <II-3> 

Equations  (II. 1)  -  (II. 3)  constitute  the  reformulated  tank  energy 
and  mass  balances  that  are  appropriate  for  volatile  cargos. 

For  a  nonvolatile  cargo  (vapor  pressure  at  normal  tempera¬ 
ture  is  much  smaller  than  atmospheric  pressure) ,  the  previous 
mass  and  energy  balances  include  many  zero  or  near-zero  terms, 
and  a  numerical  solution  using  these  equations  would  be  unduly 
complicated.  Only  the  mass  balances  (equations  (II. 2)  with  We  =  0) 
and  a  relation  that  gives  the  change  in  the  vapor  space  pressure 
when  the  vapor  space  volume  changes  are  required.  If  the  vacuum 
relief  valve  is  jammed  shut,  the  pressure  relation  follows  from 
the  adiabatic  expansion  relation: 


Ti 


-  WLi]'* 

VT  -  Vl  / 


(II. 4) 


where  the  subscript  i  indicates  the  initial  conditions  at  the 
time  the  puncture  is  created.  The  exponent  Ym  is  equal  to 
one  for  isothermal  conditions,  and  equal  to  the  ratio  of  spe¬ 
cific  heats  of  the  vapor/air  mixture  in  the  vapor  space  for 
adiabatic  conditions.  If  the  relief  valve  is  operable,  equa¬ 
tion  (II. 4)  is  still  valid  as  long  as  the  tank  pressure  re¬ 
mains  above  the  set-point  vacuum  of  the  valve.  Once  the  set- 
point  is  exceeded,  the  tank  pressure  relation  is 


FT  =  Patm  "  APvalve 


(II. 5) 


since  it  is  assumed  that  the  valve  always  remains  open  there¬ 
after.  1 


II. 2  Flow  Rate  Models 

The  next  part  of  the  venting  rate  model  relates  the  dis¬ 
charge  rates  W^,  WVQ,  and  WAq  to  the  tank  conditions. 

When  its  vapor  pressure  is  higher  than  atmospheric  pres¬ 
sure,  a  venting  liquid  may  partially  vaporize  and  a  venting  vapor 
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may  partially  condense.  These  phenomena  must  be  included  in 
the  flow  models;  the  fact  that  they  were  not  included  in  the 
models  described  in  [1]  limits  their  predictive  accuracy,  as 
will  be  discussed  further  in  Section  III. 4.  Although  there 
are  many  correlations  available  to  compute  two-phase  flow  [4] , 
they  are  all  mathematically  very  complicated  when  a  change  of 
phase  occurs  within  the  flow  stream.  For  this  reason,  a  ther¬ 
modynamic-equilibrium  type  of  flow  is  assumed,  but  an  empiri¬ 
cally  determined  factor  is  applied  to  the  effective  density 
to  increase  the  reliability  of  the  predictions. 

In  general,  a  thermodynamic- equilibrium  flow  is  given  by 

Ho  ■  CDVout  12  “-in-W  "■  2?(Jin-W,,V2 

where  the  subscript  "in"  refers  to  a  location  within  the  tank 
where  the  fluid  velocity  is  small  and  "out"  refers  to  the 
exit  plane  of  the  puncture.  The  discharge  coefficient  CD  is 
meant  to  account  for  irreversible  flow  and  other  nonideal  ef¬ 
fects.  Equation  (II. 6)  can  be  specialized  to  include  volatile 
and  nonvolatile  liquids  and  gases.  For  a  volatile  liquid,  the 
appropriate  form  is  derived  by  assuming  an  isentropic  process 
through  the  puncture;  the  result,  as  shown  in  Appendix  A,  is: 

hlo  ’  cd*lo!,lo{2  1cpl(t-V  -  V.  “  ^ 

*  (PT-PV)/PL  +  q  IZL  -  ZLh>  A 1/2  (II. 7a) 

where  o  is  the  exit-plane  density: 

liO 

°LO  "  {Vt.s  *  [aT*  <Vva-''LS»CPLln  17]/^}  <II-7b) 

Here,  the  subscript  s  indicates  the  saturated  state  at  which 
the  pressure  is  equal  to  the  pressure  outside  the  tank  (at¬ 
mospheric,  plus  the  static  head  of  water  if  the  puncture  is 
submerged) .  The  temperature  and  the  other  properties  (vLg, 
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vVg,  etc.)  can  be  computed  from  the  Chemical  Properties  File 
data.  Average  values  over  the  temperature  range  T  to  Tg 
should  be  used  for  CpL  and  o^.  The  term  (PT  -  p^)  in  equation 
(II. 7a)  allows  for  the  possibility  that  the  vapor  space  pres¬ 
sure  differs  from  the  saturation  pressure  Pv  corresponding  to 
the  liquid  temperature. 

The  factor  S  in  equation  (II. 7b)  represents  an  empirical 
way  to  incorporate  the  experimental  observation  [25]  that  the 
actual  fraction  of  the  liquid  vaporized  is  somewhat  less  than 
would  be  predicted  by  thermodynamic  equilibrium  relations. 
According  to  the  tests  described  later,  3  =  0.12  gives  a  good 
correlation  to  the  data.  It  is  recognized  that  equations 
(II. 7a)  and  (II. 7b)  are  just  an  approximation  to  a  complex 
two-phase  flow  [4] ,  but  the  comparison  of  the  model  to  the 
tests  indicates  that  the  approximation  is  satisfactory  for 
this  application. 

The  flow  relation  for  a  volatile  vapor  is  similarly: 
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(II. 8b) 


When  a  mixture  of  vapor  and  air  is  venting,  the  air  tends  to 
suppress  the  fairly  small  phase  change  effects  of  the  vapor, 
so  equation  (II. 6)  in  this  case  reduces  to  [5] : 
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and 


The  specific  volumes  of  the  vapor  and  air  in  the  tank  are 
based  on  their  respective  partial  pressures. 

For  nonvolatile  cargos,  the  flow  equations  for  vapor  or 
a  vapor/air  mixture  are  the  same  as  (II. 9a)  and  (II. 9b).  The 
flow  equation  for  discharge  of  a  nonvolatile  liquid  can  be  ob¬ 
tained  by  standard  3ernoulli-type  relations;  it  is; 

WLo  =  CDALo^L  {2[(PT-P-/PL  +  ?  (ZL-ZLh>  ]}  ^  (II‘10) 

In  all  these  models,  as  well  as  in  the  water  and  air  in¬ 
gestion  models  to  be  described  subsequently,  both  vapor  and 
liquid  outflows  can  occur  simultaneously.  This  is  handled 
analytically  by  partitioning  the  total  puncture  area  into  a 
liquid  outflow  area  ALq  and  a  vapor  outflow  area  AGq. 


Evaporation  Rate  Model 


The  discharge  rates  of  liquid  and  gas  depend  on  the  vapor 
space  pressure,  P^.  Since  evaporation  of  liquid  into  the  vapor 
space  during  the  discharge  tends  to  keep  PT  from  decreasing  as 
the  vapor  space  volume  increases,  it  is  important  to  model  the 
evaporation  realistically. 

Given  the  cargo  temperature,  the  equilibrium  pressure  of 
the  vapor  (i.e.,  the  liquid  saturation  pressure)  can  be  com¬ 
puted  from  the  Chemical  Properties  Files.  The  pressure  and 
temperature  fix  the  vapor  density,  so  equation  (II.  3)  can  be 
used  to  determine  the  evaporation  rate  required  to  maintain  the 
pressure  at  its  equilibrium  value.  AMSHAH  [1]  used  exactly  thi 
procedure.  However,  it  is  clear  that  to  create  the  required 
evaporation  rate  there  must  be  a  driving  force;  the  vapor  space 
pressure  must  be  less  than  the  liquid  saturation  pressure  [6,7] 
The  tests  conducted  during  the  present  program  confirmed  this 
conclusion.  Thus,  vy  cannot  be  related  only  to  the  cargo 


temperature,  and  equation  (II. 3)  must  be  augmented  by  an  evapora¬ 
tion  rata  model. 

Evaporation  rates  from  a  liquid  into  the  atmosphere  are 
commonly  related  to  the  difference  between  the  saturation  pres¬ 
sure  of  the  liquid  and  the  actual  partial  pressure  of  the  vapor 
in  the  air  above  the  liquid  [3] : 

We  =  3AT  (?V,sat-  V//2^T'  (II. 11) 

where  the  coefficient  a  must  be  determined  experimentally. 

For  the  present  tests,  where  evaporation  in  a  closed  container 
was  obtained,  equation  (II. 11)  could  not  be  made  to  correlate 
the  data  very  well;  further,  the  best  values  of  a  appeared  to 
vary  widely  from  test  to  test,  as  a  function  of  venting  rates 
and  chemical  properties.  By  examining  several  kinds  of  cor¬ 
relations,  it  was  found  that  a  better  data  fit  for  these  tests 
was  given  by: 

wa  *  *  1‘5  'Vsat >5//<Vl  /®’>1/4 

where  W  is  the  evaporation  rate  recuired  to  maintain  the 

e ,  sat 

vapor  partial  pressure  at  its  equilibrium  value.  The  constant 
in  equation  (11.12)  varied  from  about  1.6  to  2.2  from  test  to 
test,  with  the  best  average  being  1.9.  The  rationale  for  equa¬ 
tion  (11.12)  is  given  in  Appendix  A. 

Although  equation  (11.12)  implies  that  Pv  g  -  Pv  non" 
zero,  it  cannot  be  put  in  a  form  that  reveals  this  dependency 
explicitly.  However,  the  model  does  show  that  larger  evapora¬ 
tion  rates  give  larger  differences  between  W  -  W  ,  and 

6  ;  S21  L  6 

therefore  larger  pressure  differences,  P„  -  P„. 

In  the  computerized  solution  procedure,  the  evaporation 
rate  model  is  used  as  follows.  At  any  instant,  W  is  cal- 

6  /  Sd  w 

culated  from  equation  (II.  3)  by  assuming  that  the  vapor  is  in 
equilibrium  with  the  liquid.  Equation  (11.12)  is  then  used  to 
compute  the  corrected  Wg.  Using  this  value  of  w  ,  equation 
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(II. 3)  is  solved  for  v^,  and  finally  the  actual  vapor  pressure 
is  computed  from  the  vapor  equation  of  state:  P  =  ZPT/v  . 

V  / 

Several  iterations  are  required  for  convergence. 

When  the  pressure  at  the  puncture,  e.g.,  ?T  +  pTg(Z^-ZT^) 
is  substantially  larger  than  Pa ,  the  changes  in  ?T  during  the 
venting  have  a  relatively  small  effect  on  the  discharge.  There 
fore,  it  is  apparent  that  the  effects  of  the  evaporation  rate 
on  the  discharge  are  likely  to  be  small  for  large  tanks. 

II . 4  Air  and  Water  Ingestion  Models 

There  are  many  cases  in  which  air  or  water  ingested  into 
the  tank  through  the  puncture  can  cause  additional  cargo  to  be 
discharged.  For  example,  when  liquid  is  venting  and  the  vac¬ 
uum  relief  valve  is  jammed  shut,  a  vacuum  may  eventually  be 
created  in  the  vapor  space  that  is  sufficient  to  lower  the  tank 
pressure  at  the  puncture  to  the  ambient  pressure.  Even  though 
the  net  pressure  differential  across  the  puncture  is  now  zero, 
the  discharge  will  not  cease  if  the  liquid  level  is  still  above 
the  puncture,  although  this  was  assumed  to  be  the  case  in  [1]. 
Instead,  water  or  air  (depending  on  whether  the  puncture  is 
submerged)  will  be  ingested  to  relieve  the  vacuum.  For  sub¬ 
merged  punctures,  water  can  be  ingested  even  when  the  relief 
valve  is  operable. 

II. 4.1  Air  Ingestion  Model 

Air  can  only  be  ingested  when  (1)  liquid  is  being  vented, 
(2)  the  puncture  is  not  submerged,  and  (3)  the  static  pressure, 
PT  +  p^g  (ZL~  Z^)  ,  at  the  puncture  is  nearly  equal  to  atmo¬ 
spheric  pressure.  When  all  these  conditions  are  met,  typically 
when  the  relief  valve  is  jammed  shut,  tests  have  shown  that  the 
outflow  becomes  .intermittent.  Air  bubbles  are  ingested  to  re¬ 
lieve  the  vacuum  in  the  vapor  space ,  thereby  halting  the  out¬ 
flow  temporarily;  more  liquid  is  then  discharged,  until 
PT  +  9L5  ZLh^  a9ain  decreased  to  near  atmospheric 

pressure;  more  air  bubbles  are  ingested,  and  so  on,  until  the 
liquid  level  falls  below  the  top  of  the  puncture;  see  Figure 
II. 2a. 
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FIGURE  II. 2.  AER  AND  WATER  INGESTION  MODELS 


The  available  literature  on  this  kind  of  intermittent 
flow  is  sparse  [9] .  A  relatively  simple  but  plausible  model 
has  therefore  been  developed  for  the  venting  rate  model.  The 
model  contains  one  empirical  parameter  (the  bubble  size  fac¬ 
tor)  that  is  selected  to  give  agreement  with  test  results. 

In  summary,  the  model  is  described  as  follows: 

o  Air  ingestion  begins  when  the  combination  of  tank 

liquid  level  and  vapor  space  pressure  has  decreased 
nearly  to  atmospheric  pressure: 

PT  ‘  Patm  +  *pL  <ZL  "  *ih>  1  °-001Patm 

o  When  (11.13)  is  satisfied,  an  air  bubble  is  in¬ 
gested  whose  volume  is  a  specified  multiple  (the 
bubble  size  factor)  of  a  sphere  having  a  cross- 
sectional  area  equal  to  the  puncture  area, 
o  The  vapor  space  pressure  is  increased  by  an  amount 
4Pt  corresponding  to  the  addition  of  this  bubble  to 
the  vapor  space;  for  a  volatile  cargo,  the  APT  cal¬ 
culation  is  accomplished  by  integrating  equation 
(II. 1)  over  an  infinitesimal  time  step,  At,  with 
W^At  set  equal  to  the  ingested  bubble  mass;  for  a 
nonvolatile  cargo,  the  change  in  the  specific  vol¬ 
ume  of  the  vapor  space  gas  mixture  is  computed,  and 
APt  is  then  calculated  by  ideal  gas  laws, 
o  Since  the  left  hand  side  of  equation  (11.13)  is  now 
larger  than  the  right  hand  side,  additional  liquid 
is  discharged,  in  accordance  with  the  appropriate 
energy  balances  and  flow  relations;  the  outflow 
continues  until  equation  (11.13)  is  again  satisfied. 
Another  air  bubble  is  ingested,  and  the  process  is 
repeated. 


o 
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The  right  hand  side  of  equation  (11.13)  is  set  equal  to 
0.001  Patm,  rather  than  to  zero,  to  avoid  numerical  problems 
in  the  computerized  model,  by  ensuring  that  the  outflow  has 
not  ceased  completely  before  the  bubbles  are  ingested.  Experi¬ 
ments  have  indicated  that  the  value  of  the  right  hand  side  is 
between  zero  and  about  0.002  ?atm>  so  equation  (11.13)  is  rea¬ 
sonably  accurate.  As  the  vapor  space  volume  increases,  the 
increase  in  the  pressure  AP^  due  to  the  ingestion  of  a  single 
bubble  becomes  smaller,  and  the  net  AP  driving  the  liquid  out¬ 
flow  therefore  becomes  smaller;  further,  the  hydrostatic  head 
ZL  “  ZLh  sma^er  than  it  was  originally.  To  prevent  an  ex¬ 
cessive  number  of  computational  iterations,  more  than  one  bubble 
can  be  admitted  when  necessary.  The  criterion  used  in  the  com¬ 
puterized  model  is  that  the  value  of  PT  +  p^g  (Z^-  Z^h)  after 
each  bubble  (or  bubbles)  has  been  admitted  should  be  reasonably 
close  to  the  value  obtained  after  the  very  first  bubble  is  in¬ 
gested.  This  procedure  also  tends  to  bring  the  vapor  space 
pressure  smoothly  to  atmospheric  pressure  just  as  the  puncture 
is  uncovered. 

By  comparison  of  the  venting  rate  model  to  the  air  inges¬ 
tion  tests,  the  best  value  of  the  bubble  volume  factor  mentioned 
above  has  been  found  to  be  5.6.  This  gives  an  effective  bubble 
diameter  of  about  1.8  times  the  puncture  diameter,  which  agrees 
qualitatively  with  visual  observations. 

For  very  large  punctures ,  and  especially  for  long  vertical 
punctures,  the  outflow  may  not  completely  cease  during  the  air 
ingestion;  that  is,  the  bubbles  may  be  ingested  only  over  the 
top  part  of  the  puncture.  This  kind  of  ingestion  was  not  ob¬ 
served  in  the  smaller  punctures  used  in  the  present  tests,  but 
if  it  does  occur,  the  model  described  herein  will  provide  a 
first  approximation  to  the  actual  outflow  rate. 

II. 4. 2  Water  Ingestion  Models 

When  the  puncture  is  below  the  water  line,  the  tests  have 
shown  that  water  will  be  ingested  whenever  the  net  pressure 
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difference  across  the  puncture  approaches  zero.  Once  started, 
water  ingestion  continued  until  the  puncture  was  under  water 
on  both  sides.  In  contrast  to  air  ingestion  where  the  outflow 
became  intermittent,  here  the  outflow  remained  steady  but 
shared  the  puncture  area  with  the  water  inflow.  In  other 
words,  there  was  a  water  stream  entering  the  tank  and  simul¬ 
taneously  a  cargo  liquid  stream  leaving  the  tank. 

The  models  developed  here  assume  the  cargo  and  the  in¬ 
gested  water  do  not  mix  or  react  chemically.  This  is  a  con¬ 
servative  assumption  for  many  cargos,  since  the  venting  rate 
will  be  predicted  to  be  larger  than  will  actually  occur.  In 
any  event,  a  model  that  included  mixing  or  chemical  reactions 
in  the  tank  would  be  excessively  complicated.  Since  the  cargo 
and  the  liquid  do  not  mix,  the  ingested  water  will  displace 
a  lighter  cargo  upward,  thereby  exposing  more  cargo  to  the 
puncture  and  increasing  the  static  pressure  at  the  puncture; 
see  Figure  II. 2b.  When  the  water  is  lighter  than  the  cargo, 
the  ingested  water  will  float  on  top  and  thereby  increase  the 
static  pressure  at  the  puncture;  see  Figure  II. 2c.  In  both 
cases,  the  ingested  water  allows  additional  cargo  to  be  dis¬ 
charged. 

The  driving  force  for  water  ingestion  is  the  buoyancy 
(negative  or  positive)  of  the  water  relative  to  the  cargo, 
which  forces  the  water  into  the  tank  when  the  cargo  outflow 
momentum  is  low  (i.e.,  when  the  net  positive  pressure  differ¬ 
ence  across  the  puncture  is  near  zero) .  The  tests  showed  that 
the  water  entering  the  tank  under  these  conditions  soon  broke 
up  into  a  stream  of  discrete  droplets.  Since  an  analysis  of 
water  ingestion  could  not  be  found  in  the  literature,  several 
r.ew  but  approximate  models  were  developed.  The  models  contain 
one  empirical  constant  that  can  be  chosen  to  give  better  corre¬ 
lation  to  the  test  data.  In  summary,  they  are  described  as: 
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o  Tank  conditions  for  which  the  pressure  differential 
across  the  puncture  is  zero  are  determined  by  an 
equation  similar  to  equation  (11.13)  of  the  air  in¬ 
gestion  model,  but  with  the  right  hand  side  replaced 
by  the  hydrostatic  head  gpw  (Zw~  ZLh)  . 

o  Ingestion  begins  just  prior  to  this  condition  of 

zero  net  A?,  such  that  a  small  positive  pressure 
differential,  the  magnitude  of  which  is  determined 
with  the  aid  of  the  empirical  constant  mentioned 
earlier,  will  continue  to  force  out  cargo. 

o  The  water  inflow  remains  in  a  stream  tube;  this 
assumption  allows  the  calculation  of  a  AP  acting 
on  the  water  in  the  direction  to  force  water  into 
the  tank. 

o  The  puncture  area  is  shared  between  the  cargo  out¬ 
flow  and  the  water  inflow  such  that  the  inflow  and 
outflow  rates  maintain  the  same  pressure  difference 
driving  the  cargo  outflow  as  at  the  start  of  inges¬ 
tion. 

When  the  water  is  lighter  than  the  cargo  and  the  cargo 
static  pressure  at  the  puncture  is  equal  to  the  outside  pres¬ 
sure  (atmospheric  plus  water  hydrostatic  head) ,  the  net  pres¬ 
sure  difference  from  the  "inlet"  to  the  "outlet"  of  the  water 
streamtube  mentioned  above  (Figure  II. 2c)  can  be  calculated  by 
considering  the  various  static  heads;  the  result  is; 

APW  *  9  (pL-PW}  (11-14) 

which  is  just  the  buoyancy- induced  pressure  on  a  vertical 
water  column  of  equal  length.  When  the  water  is  heavier  than 
the  cargo,  the  pressure  differential  is  similarly; 

Apw  "  g  (pw-ol}  (zLh'zwr)  (ii.i5) 

In  both  cases,  the  water  inflow  rates  are  calculated  by: 


(11.16) 


W. 


w 


CDAinpW  ' 2  “PW/pW 


where  A^n  is  the  part  of  the  puncture  occupied  by  the  water 
stream.  There  is  little  or  no  change  of  phase  (evaporation) 
of  the  liquid  outflow  during  water  ingestion  since  the  tank 
pressure  is  nearly  in  equilibrium  with  the  outside  pressure; 
see  equations  (II. 7b)  and  (II. 8b).  Thus,  for  simplicity,  the 
cargo  venting  rate  is  calculated  using  an  incompressible  flow 
relation.  The  water  inflow  rate  is  supposed  to  be  just  large 
enough  to  maintain  the  original  net  AP  at  the  time  when  inges¬ 
tion  started,  so  the  cargo  outflow  rate  is: 

_  -sl/2 

W  =  C  A  o  [2  RPm-  -  P-  _  )/p_  +  g  (z_.  -  ZT  _  )  /K>  (11.17) 

Lo  U  OUtpL  \  |_ ‘  Tl  T,Eq./ML  *  Li  L.Eq.J/  J 

where  P_  _  and  Z  are  the  values  of  P_  and  Z_  which  first 

1  t  •  Lt  f  £G [  •  1  Xj 

satisfy  the  water  ingestion  form  of  equation  (11.13).  PTi  and 
are  values  at  some  slightly  earlier  time;  in  the  computer¬ 
ized  model,  these  correspond  to  the  values  at  the  time  step  pre¬ 
ceding  the  start  of  water  ingestion.  The  actual  &P  across  the 
puncture  may  be  only  a  fraction  of  the  value  (pTj_  ~  PT  Eg  )/pL 
+  g  ( Z_ .  -  Z_  )  selected  by  the  computer  scheme.  Thus,  an 
empirical  constant  K  is  inserted  in  equation  (11.17);  the  best 
value  obtained  from  the  test  comparisons  is  K  =  83.  It  will  be 
shown  later  that  there  is  reason  to  believe  that  K  =  83  is  a 
universal  value,  and  not  just  a  value  specific  to  a  given  test 
and  computer  prediction. 

When  the  ingested  water  is  lighter  than  the  cargo,  it  floats 
on  the  liquid  surface  and  suppresses  any  further  evaporation. 

When  it  is  heavier  and  thus  sinks,  evaporation  is  possible,  but 
it  is  neglected  in  the  models  because  the  energy  balances  needed 
to  predict  temperature  changes  become  much  more  complicated. 

In  effect,  the  cargo  is  treated  as  an  incompressible  liquid  and 
a  non-condensable  gas  after  water  is  ingested. 

As  mentioned  earlier,  the  inflow  rate  of  water  is  adjusted 
so  as  to  maintain  the  pressure  differential  driving  the  cargo 
outflow.  When  the  water  is  heavier  than  the  cargo  (Figure 
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II. 2b),  this  merely  requires  that  the  two  volumetric  flow  rates 
be  equal  so  that  ZT  remains  constant;  in  other  words 


W 


w 


=  p  WT  /p. 
w  Lo'  L 


(11.18) 


Equations  (11.16),  (11.17),  and  (11.18)  are  all  made  compatible 
by  allocating  the  puncture  area  between  the  two  streams  as; 

A.  +  A  .  =  A  (11.19) 

m  out  p  '  J 

These  equations  are  sufficient  to  determine  all  the  flow  rates. 


When  the  water  is  lighter  than  the  cargo  (Figure  II. 2c) , 
the  relation  between  W  and  w_  depends  on  whether  the  relief 
valve  is  operable.  When  it  is,  the  vapor  space  pressure  remains 
constant.  Thus,  to  maintain  a  constant  a?  across  the  puncture, 
the  weight  of  the  water  ingested  must  equal  the  weight  of  the 
cargo  discharged,  for  otherwise  the  hydrostatic  head  above  the 
puncture  would  change.  The  flow  relation,  consequently,  is 


w. 


w 


=  w. 


Lo 


(II. 20a) 


On  the  other  hand,  when  the  relief  valve  is  jammed,  the  condi¬ 
tions  specified  by  equation  (II. 20a)  would  compress  the  vapor 
space  gas  since  a  larger  volume  of  water  would  be  admitted  than 
the  volume  of  cargo  discharged.  Thus,  the  increase  in  PT  and 
the  hydrostatic  pressure  changes  must  be  balanced  in  order  to 
hold  constant  the  AP  driving  the  water  outflow.  The  appropri¬ 
ate  relation  is  derived  by  expressing  the  AP  at  the  puncture 
in  terms  of  the  vapor  space  compression  and  the  changes  in 
water  and  cargo  levels;  the  result,  as  shown  in  Appendix  A,  is: 


W. 


w 


W. 


Lo 


g  (zT  -  zm)  +  rmpT//pL 

g  (ZT~ZWT)  +  rmPT,/pW 

*/ 


(II. 20b) 


Equations  (11.13)  -  (11.20)  constitute  the  models  of  water 
ingestion  for  those  cases  where  some  cargo  liquid  is  discharged 
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prior  to  the  start  of  ingestion.  There  are  also  situations 
where  water  can  be  ingested  immediately  upon  puncturing  the 
tank.  This  would  occur  whenever 


PTi  "  Ptt»  *  «U-2Lh>  '  (JW  ‘  *lh>  (11.21) 

For  these  conditions,  cargo  would  never  be  discharged  unless 
water  ingestion  permits  it. 

When  ow  <  s^,  it  is  assumed  that  water  will  enter  the 
tank  and  rise  to  the  top  without  any  cargo  being  discharged, 
until  the  tank  pressure  at  the  puncture  is  equal  to  the  out¬ 
side  pressure: 

PT  “  Patm  +  gpW  (ZWT,£q.  ~  ZLi)  +  g3L  (ZLi  “  ZLh5 

“  gpw  (Zw~  2^)  (11.22) 

Of  course,  if  the  required  level  of  water  is  greater 

than  the  elevation  of  the  top  of  the  tank,  the  tank  will  fill 
up  with  water  before  equation  (11.22)  is  satisfied;  this  case 
is  modeled  separately.  If  Z„_  _  <  Z„,  water  will  continue 

Wi  f  •  A, 

,to  flow  into  the  tank,  but  now  cargo  will  be  discharged  since 
there  is  a  net  positive  dP  across  the  puncture.  At  this  point, 
the  venting  model  is  the  same  as  for  the  similar  water  inges¬ 
tion  models  described  earlier,  although  the  cargo  outflow  rate 
equation  can  be  simplified,  as  shown  in  Appendix  A,  to  give: 

WLo  -  CDA0ut»i  (IX. 23) 

ZWT,Eq.  can  exceed  zx  when  tile  valve  is  operable, 

and  does  not  increase  as  the  water  enters.  When  this  occurs, 
it  is  assumed  that  the  water  completely  fills  up  the  vapor 
space.  The  pressure  thereafter  in  the  vanishingly  small  vol¬ 
ume  between  the  top  of  the  tank  and  the  water  is  assumed  to 
remain  constant  at  the  initial  value  maintained  by  the  relief 
valve.  Thus,  as  further  water  enters  in  an  attempt  to  satisfy 


21 


equation  (11.22),  cargo  will  be  forced  out  at  the  sane  volumet¬ 
ric  flow  rate.  For  this  case,  there  is  a  net  AP  across  the 
puncture  forcing  water  into  the  tank,  so 


VI 


w 


and 


CDAinpW  ^  (Patm  "  PTi  ^  PW  +  9  (ZW  “  Zt'! 
-  g  (pl-pw)  <Z-L-ZLh )l}1/2 
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(II. 24a) 


(II. 24b) 


There  are  no  other  relations  available  to  determine  the  way  the 

Duncture  area  is  shared,  so  it  is  assumed  that  A.  *  A  =A  /2. 

m  out  p 

When  the  relief  valve  is  jammed  shut,  the  vapor  space 
pressure  will  increase  as  the  water  enters,  and  the  condition 
that  >_  ZT  is  impossible.  For  this  case,  equation 

(11.22)  must  be  solved  simultaneously  with  the  relation  be¬ 
tween  vapor  space  volume  and  pressure: 


(11.25) 


where  the  subscript  "i"  denotes  the  initial  conditions. 

When  ow  >  pl,  the  water  initially  ingested  into  the  tank 
will  lift  the  cargo  until  there  is  a  hydrostatic  balance  at 
the  puncture.  Since  the  original  depth  of  cargo  is  ZL^,  the 
equilibrium  relation  is 


PT  "  Patm  *  gpL  (ZLi  +  ZWT,Eq.  ”  ZLh)  *  gp  W  ( 2W  “  ZLh )( 11  *  26  } 

where  ZLi  +  ZWT  ^  is  the  new  elevation  of  the  cargo  surface. 
Assuming  for  the  monvent  that  (1)  Z_  .  +  Z._  _  <  Z_  (i.e.,  the 

IjX  WT  f  •  X 

new  liquid  level  is  below  the  top  of  the-  tank)  ,  and  (2)  the 
water  in  the  tank  does  not  cover  the  vent,  further  water  will 
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be  ingested,  but  now  cargo  will  be  discharged  at  the  same 

volumetric  rate.  The  model  for  this  kind  of  venting  is  the 

same  as  the  analogous  models  described  earlier,  except  again 

it  must  be  assumed  that  A.  =A  .  =  A  /2.  Ingestion  will  con- 

in  out  p  3 

tinue  until  the  puncture  is  covered  by  water  on  both  sides. 

If  the  water  initially  covers  the  vent,  cargo  is  never  dis¬ 
charged.  If  zLi  +  ZWT  Eq  >  zt'  w^^-ch  can  occur  only  when  the 
relief  valve  is  operable,  it  is  assumed  that  the  cargo  is 
lifted  entirely  to  the  top  of  the  tank  and  the  vapor  space 
pressure  remains  at  its  initial  value.  Then,  as  further  water 
is  ingested,  cargo  is  discharged  at  the  same  volumetric  flow 
rate.  The  water  inflow  rate  is  then  given  by 

w„  -  CDAin“w{2  +  9  (zw-2Lh> 

-  g  („L/o„)  (ZT-2Lh))}1/2  (11.27) 

which  is  based  on  the  net  AP  from  inside  to  outside  at  the 
center  of  the  puncture.  When  the  relief  valve  is  jammed  shut, 
an  expression  analogous  to  equation  (11.25)  is  used  to  compute 

ZWT , Eq .  ' 

This  completes  all  the  water  ingestion  models. 

II . 5  Choked  Flow  Models 

Some  cargos  (e.g.,  butane)  have  vapor  pressures  at  normal 
temperatures  that  are  substantially  higher  than  atmospheric 
and  are  carried  in  pressurized  tanks.  If  one  of  these  tanks 
is  punctured,  the  pressure  differential  across  the  opening  may 
be  so  large  that  the  outflow  velocity  is  limited  by  the  speed 
of  sound  [5];  in  other  words,  the  flow  is  "choked"  to  some 
lower  value  than  that  predicted  by  the  models  described  in 
II. 2,  all  of  which  assume  that  the  pressure  in  the  outflow 
stream  at  the  puncture  exit  is  the  outside  ambient  pressure. 
For  these  cases,  choked  flow  models  are  developed. 

For  venting  of  a  gas,  the  limiting  exit  pressure  for 
choked  flow  is  (5] : 
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If  the  calculated  Pc  is  greater  than  ambient  (atmospheric, 
plus  water  hydrostatic  head  if  the  puncture  is  submerged) ,  the 
outflow  will  be  choked,  and  P  must  be  used  in  equation  (II. 9a) 

v 

instead  of  Pb  to  compute  WVq.  As  the  vapor  continues  to  vent, 
PT  will  gradually  decrease  until  eventually  Pc  =  Pao;  at  this 
point,  the  outflow  is  no  longer  choked,  and  equation  (II. 9a) 
can  be  used  directly.  Note  that  equation  (11.28)  considers 
the  effects  of  two-phase  flow  for  vapor  as  being  negligible. 

Ordinarily,  flow  of  a  liquid  cannot  be  choked  because  of 
the  very  high  speed  of  sound  for  liquids.  For  a  volatile 
liquid,  however,  the  change  in  phase  of  the  outflow  can  easily 

lower  the  speed  of  sound  to  a  fairly  small  value.  The  exit 
conditions  for  choked  flow  of  a  volatile  liquid  can  be  deter¬ 
mined  by  maximizing  the  mass  flow  per  unit  area: 

• 

ClX"  =  [  2  (hL  “  hLc  "  XcXc)  1  ^ /[vLc  +  Sxc  (VVc  "  VLc)]  (11.29) 

D  O 


where  xc  represents  the  relative  mass  fraction  of  vapor  in 
the  outflow.  (The  static  head  above  the  puncture,  pgT  (ZT  -  ZT.  )  , 
is  neglected  in  this  equation,  since  it  should  be  small  com¬ 
pared  to  the  tank  pressure  Pv  whenever  choked  flow  is  likely.) 
The  maximum  of  equation  (11.29)  corresponds  to 


hL~hLc~Vc 

VLc  +  8XC  (VVC-VLC> 


l 


3TT  (h 


Lc  +  xcV 


ctT  tv 


Lc  +  SXc  <vV=-vLc» 


(11.30) 


Equation  (11.30)  can  be  solved  numerically  using  the  known 
cargo  pressure  and  temperature  to  determine  the  choked  flow 
exit  temperature  T  and  corresponding  saturation  pressure,  P  . 

w  w 

If  this  pressure  is  greater  than  Pa,  the  flow  is  choked,  and 
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these  values  of  Tc  and  Pc  must  be  used  in  equations  (II. 7)  to 
compute  WLq.  As  the  liquid  vents,  PT  may  eventually  decrease 
to  the  point  where  the  exit  pressure  Pc  required  to  choke  the 
flow  is  smaller  than  Pa;  after  this  point,  equations  (II. 7) 
can  be  used  directly.  Depending  on  the  saturation  pressure  of 
the  liquid,  however,  it  is  entirely  possible  for  the  liquid 
flow  to  remain  choked  throughout. 

The  choked  flow  model  given  by  equations  (11.29)  and 
(11.30),  although  similar  to  that  described  in  the  Interim 
Report  [3] ,  assumes  equilibrium  flow  rather  than  attempting 
to  incorporate  nonequilibrium  effects  [10,11]  as  suggested  in 
[3].  The  assumption  of  equilibrium  flow  is  in  general  agree¬ 
ment  with  the  assumptions  used  for  nonchoked  flows. 

In  the  computer  code,  equation  (11.30)  is  solved  by  ex¬ 
panding  the  left  hand  side  in  powers  of  AT,  where  Tc  =  T-AT. 

The  derivatives  on  the  right  side  are  evaluated  numerically 
by  assuming  two  values  for  Tc.  As  shown  in  Appendix  A,  the 
temperature  difference  AT  is  then  found  from: 


where  [D]  is  the  computed  approximate  value  of  the  derivatives: 


(II. 31b) 


The  parameter  E  is  related  to  the  saturation  pressure  correla¬ 
tion  (see  equation  (11.33)  of  Section  II. 6)  and  is: 


E  =  [BqT/  (CQ  +  T* )  2  ]  in  10  -  1.5 


(II.31c) 
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(T*  is  in  °C,  to  conform  to  the  Chemical  Properties  File  usage; 
all  the  other  temperatures  are  °K.)  One  or  two  iterations  of 
equations  (11.31)  have  been  shown  to  be  sufficient  to  deter¬ 
mine  T  . 
c 

II . 6  Review  of  Chemical  Properties  File 

The  previous  technical  review  [2]  indicated  that  the  ther¬ 
modynamic  data  in  the  Chemical  Properties  File  were  in  error 
in  some  cases  and  the  temperature  range  covered  by  the  corre¬ 
lating  equations  was  not  sufficient  in  other  cases.  These 
conclusions  have  been  confirmed  by  the  present  review.  Six 
representative  chemicals,  all  of  which  might  be  transported 
in  a  state  where  phase  changes  can  occur  so  that  thermodynamic 
properties  are  important,  were  reviewed  in  detail  to  document 

any  errors;  these  chemicals  are:  anhydrous  ammonia,  butane, 
chlorine,  LNG,  methyl  alcohol,  and  propane. 

Based  on  this  review  of  the  representative  chemicals,  the 
following  conclusions  have  been  drawn  about  the  adequacy  of 
the  present  Chemical  Properties  File  for  use  in  the  venting 
rate  model. 

1.  Correlating  equations  for  liquid  density  and  spe¬ 
cific  heat  are  slightly  inaccurate,  but  probably 
acceptable . 

2.  Correlating  equations  for  saturated  vapor  proper¬ 
ties  are  probably  not  acceptable;  in  particular, 
the  specific  heat  for  saturated  vapor  (not  super¬ 
heated  vapor)  and  the  compressibility  factor  should 
be  included. 

3.  Correlating  equations  for  saturation  pressure  as  a 
function  of  temperature  are  moderately  inaccurate 
and  can  lead  to  inaccurate  predictions  of  tank 
pressure  during  venting. 

4.  The  single  value  of  latent  heat  given  for  each 
chemical  differs  by  as  much  as  20%  from  the  true 
value  over  the  range  of  applicable  temperatures 
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needed  in  the  venting  rate  model  energy  balance 
and  flow  equations. 

Recommended  correlations  for  the  six  reviewed  chemicals, 
based  on  the  indicated  data  sources,  are  given  in  Appendix  B. 

The  recommended  correlations  for  latent  heat  are  derived  from 
the  Clausius-Clapeyron  equation,  X  =  (T*+  273.17)  (vv  -  v.)  (dPv/dT) , 
using  the  new  correlations  given  for  vapor  compressibility, 
liquid  density,  vapor  compressibility  factor,  and  vapor  pres¬ 
sure.  Thus,  a  separate  correlation  for  X  is  not  needed,  since 
it  can  be  computed  from  the  included  correlations  by  the  rela¬ 
tion 


X  =  ZRBq  Zn  10  [ (T*  +  273.2)/ (T*  + CQ) 3 2 

-  5.501  x  10~4  B  Ptr  (T*  +  273.2)  /  (T*  +  C. )  p  T  (11.32) 

O  V  O  Li 

Here,  Bq  and  CQ  are  constants  in  the  vapor  pressure  correlation: 

i°g10pv  -  Ao  '  V(T‘  +  co>  (II-33) 

(In  both  equations  (11.32)  and  (11.33),  the  units  of  tempera¬ 
ture  T*  are  °C.) 

Nonvolatile  cargos  are  usually  transported  at  near  atmo¬ 
spheric  pressures.  Venting  of  gas  is  therefore  unlikely  to  be 
important,  and  there  will  be  little  or  no  change  of  phase  dur¬ 
ing  liquid  venting.  Under  these  restrictions,  the  liquid  den¬ 
sity  is  the  only  cargo  property  that  is  important.  The  present 
Chemical  Property  File  data,  consequently,  are  sufficiently 
accurate  for  these  cases. 


III.  EXPERIMENTAL  PROGRAM 


There  were  three  primary  objectives  of  the  test  program: 

1.  Obtain  data  to  evaluate  the  assumptions  of  the 
analytical  models 

2.  Gain  insight  into  air  and  water  ingestion 

3.  Obtain  discharge  rate  data  for  a  variety  of  test 
conditions  to  validate  the  overall  reliability  of 
the  venting  rate  model. 

III.l  Scaling  Considerations 

All  the  tests  were  conducted  using  tanks  that  were  much 
smaller  than  full  scale.  It  was  important,  then,  that  the 
simulated  cargos  and  the  initial  tank  conditions  be  selected 
so  as  to  fairly  represent  full  scale  venting  behavior.  The  liq¬ 
uid/vapor  properties  and  the  test  conditions  needed  to  achieve 
the  desired  simulation  were  determined  from  a  similitude  analy¬ 
sis  [12]  in  which  the  important  cargo  and  tank  parameters  were 
arranged  in  various  nondimens ional  groups.  This  analysis  showed 
clearly  that  an  exact  simulation  would  be  difficult  to  obtain. 
For  example,  the  static  head  of  liquid  in  the  tank,  gpTZ_,  was 
reduced  considerably  in  the  tests  in  comparison  to  full  scale. 

To  obtain  similarity  between  the  pressure  differentials  driving 
the  discharges  would  have  required  a  reduction  of  both  the  tank 
vapor-space  pressure  and  the  atmospheric  pressure  in  direct  pro¬ 
portion  to  goLZL.  Although  it  would  have  been  feasible  to  use 
a  simulated  cargo  whose  vapor  pressure  was  less  than  full  scale 
in  accordance  with  this  scaling  criterion,  conducting  the  entire 
test  in  a  vacuum  chamber  or  other  apparatus  to  reduce  the  atmo¬ 
spheric  pressure  would  have  been  extremely  difficult  and  expen¬ 
sive.  For  this  reason,  tests  achieving  exact  similarity  to  some 
given  prototype  were  not  attempted.  The  correct  overall  phe¬ 
nomena  were,  however,  represented  fairly  well,  as  will  be  de¬ 
scribed  in  the  next  section.  Lack  of  exact  similarity  does  not 
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seriously  limit  the  validation  of  the  analytical  model  because 
the  model  can  be  used  to  predict  results  for  the  conditions 
used  in  the  tests  and  the  predictions  compared  with  these  test 
results . 

III. 2  Summary  of  Tests 

III. 2.1  Test  Series  No.  1:  Discharge  Coefficients 

The  original  venting  rate  model  [1]  incorporated  a  fixed 
value  of  0.8  for  the  discharge  coefficient  of  the  puncture. 
Since  the  predicted  venting  rate  is  directly  proportional  to 
the  discharge  coefficient,  the  first  test  series  was  designed 
to  evaluate  the  implicit  assumption  that  there  is  a  negligible 
variation  of  CD  with  puncture  geometry,  cargo  properties,  and 
discharge  velocity. 

The  tests  were  conducted  with  an  open  cylindrical  steel 
tank  (Figure  III. la)  having  a  diameter  of  57.2  cm  (22.5  in.), 
height  of  87.6  cm  (34.5  in.),  and  total  capacity  of  224.8  li¬ 
ters  (59.4  gallons).  A  circular  hole  approximately  61  cm  in 
diameter  was  machined  in  the  wall  of  the  tank  near  the  bottom, 
around  which  a  "picture  frame"  fixture  was  mounted  flush  to 
the  outer  surface.  Simulated  punctures  were  cut  in  sheet 
metal  of  about  the  same  thickness  as  the  tank  wall  (0.16  cm 
or  1/16  in.)  and  attached  to  the  fixture,  as  indicated  in  Fig¬ 
ure  III. la.  In  this  way,  three  different  puncture  shapes 
could  be  readily  tested  using  the  same  setup: 

o  circular 

o  horizontal  rectangular  (width/height  =  5) 

o  vertical  rectangular  (height/width  =  5) 

Each  puncture  shape  was  tested  with  four  different  edge  condi¬ 
tions  : 

o  smooth 

o  jagged  petals,  all  pushed  inward  uniformly 

o  jagged  petals,  all  pushed  outward  uniformly 

jagged  petals,  pushed  inward  randomly 


o 


87.6  cm 


LIQUID 


PUNCTURE  MOUNTING 
FRAME  AND  "FLY 
AWAY"  DOOR 


b)  TYPICAL  SIMULATED  PUNCTURE  WITH  JAGGED  EDGE 

FIGURE  III.l.  TEST  APPARATUS  FOR  DETERMINATION  OF 
DISCHARGE  COEFFICIENTS 
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The  smooth-edge  punctures  served  as  a  baseline  with  which  to 

compare  the  more  realistic  jagged  punctures.  A  photograph  of 

a  typical  jagged  edge  puncture  is  shown  in  Figure  III. lb. 

Except  for  some  preliminary  tests,  the  open  area  of  all  the 

2 

punctures  was  11.4  cm  ,  corresponding  to  circular  holes  of 
diameter  3.8  cm  (1.5  in.)  and  to  rectangular  holes  7.55  cm  by 
1.51  cm  (2.97  in.  by  0.59  in.).  To  vary  the  Reynolds  number 
of  the  discharge  over  a  wide  range,  three  different  nonvola¬ 
tile  liquids  were  used:  water,  aqueous  glycerine  solution, 
and  hydraulic  oil.  The  densities  and  viscosities  of  these 
liquids  are: 


LIQUID 


DENSITY  (g/cm3) 


VISCOSITY  (cp) 


Water  0.998 
40%  Water  -  60%  Glycerine  1.157 
Hydraulic  Oil  0.878 


0.98 

8.63 

64.90 


For  a  given  test,  the  tank  was  filled  to  a  preselected 
level,  generally  64.8  cm  (35.5  in.)  above  the  centerline  of 
the  puncture.  The  discharge  was  started  by  opening  a  quick- 
release  valve.  Liquid  levels  in  the  tank  were  measured  visu¬ 
ally  during  the  test  by  a  graduated  sight  tube  attached  to  the 
tank,  and  flow  times  were  measured  with  a  stop  vTatch.  The 
tests  were  terminated  when  the  liquid  level  was  still  above 
the  puncture  (generally  14.6  cm,  or  5.75  in.).  Each  test  was 
repeated  at  least  three  times.  Test-to-test  correlation  was 
excellent,  with  the  total  discharge  time  never  varying  by  more 
than  0.5  second. 

As  is  the  usual  practice  [13]  ,  the  discharge  coefficient 
was  calculated  by  using  Torricelli's  equation  to  derive  a  rela¬ 
tion  between  CD  and  the  test  measurements: 


C 


D 


/2H“7g 


(III.l) 


where  and  are  the  initial  and  final  liquid  depths  in  the 
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tank,  and  Lt  is  the  total  discharge  time.  A  representative 
Reynolds  number  for  the  discharge  is  given  by  p.  /g  (H.  +  H-) A '/y  , 
which  is  based  on  the  discharge  velocity  corresponding  to  the 
average  liquid  depth. 

Summaries  of  the  test  results  are  given  in  Table  III.l. 

The  discharge  coefficient  varied  much  more  with  geometry  than 
with  Reynolds  number,  at  least  for  the  Reynolds  numbers  that 
are  of  most  interest.  In  particular,  the  horizontal  or  rec¬ 
tangular  punctures  with  outward-pushed  petals  give  CD's  that 
are  25%  larger  than  for  the  other  puncture  geometries.  The 
test  values  of  CD  for  the  smooth-edge  circular  puncture  agree 
reasonably  well  with  published  data  [13]  which  show  that  CD 
should  decrease  from  about  0.68  to  about  0.61  as  the  Reynolds 
number  increases  from  about  1000  to  over  100,000  or  more;  this 
favorable  comparison  validates  the  experimental  procedure. 

To  summarize,  these  data  show  that  a  discharge  coeffi¬ 
cient  of  0.65  should  be  adequate  in  practice  for  most  punc¬ 
tures.  The  single  exception  to  this  recommendation  would  be 
a  rectangular  puncture  having  ragged  edges  pointing  predomi¬ 
nately  outward;  CD  for  such  a  puncture  averages  about  0.825 
for  large  Reynolds  numbers.  Note  that  the  value  of  CD  *  0.8 
used  in  the  AMSHAH  model  [1]  roughly  corresponds  to  the  singu¬ 
lar  exception  of  the  tests. 

The  wall  thicknesses  of  the  tanks  used  in  Test  Series  2 
and  3  were  comparable  to  the  hole  diameters;  for  these  holes, 

CD  was  found  to  be  about  0.73.  It  is  believed  that  these  kinds 
of  punctures  are  not  relevant  to  actual  ship  tank  punctures, 
and  so  the  data  are  not  included  in  Table  III.l. 

III. 2. 2  Test  Series  No.  2  :  Air  and  Water  Ingestion 
Studies" 

The  second  series  of  tests  was  designed  both  to  obtain 
insight  about  liquid  discharge  during  air  or  water  ingestion 
and  to  obtain  quantitative  data.  In  order  to  make  visual  ob¬ 
servations,  the  tank  was  constructed  of  transparent  acrylic 
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plastic.  The  tank,  as  sketched  in  Figure  III. 2,  was  a  closed 
cylinder  27.6  cm  (10.87  in.)  in  diameter  and  61.0  cm  (24  in.) 
high,  with  a  total  capacity  of  36.5  liters  (9.65  gallons). 

Two  simulated  circular  punctures  were  drilled  into  the  tank 
wall  4.4  cm  (1.75  in.)  above  the  bottom.  The  smaller  puncture 
was  2.2  cm  (0.875  in.)  in  diameter,  and  the  larger  was  3.2  cm 
(1.25  in.).  Discharge  coefficients  for  these  punctures  were 
determined  to  be  0.72  -  0.73,  values  higher  than  might  be  ex¬ 
pected  from  Table  III.l,  but  which  are  the  result  of  the  much 
larger  wall  thickness  (about  1.2  cm)  of  this  tank. 

Air  Ingestion  Studies.  For  these  tests,  water  was  used 
as  a  convenient  simulation  of  a  nonvolatile  cargo.  Several 
preliminary  tests  were  conducted  to  determine  what  vapor  space 
vacuum  pressures  during  the  discharge  could  be  held  nearly 
constant  by  manually  adjusting  the  simulated  relief  valve  (Fig¬ 
ure  III. 2)  . 

The  initial  set  of  data-tests  was  conducted  to  determine 
the  functional  dependency  of  the  starting  point  of  air  inges¬ 
tion  with  liquid  level  and  vapor  space  pressure.  Within  the 
accuracy  of  the  measurements  of  liquid  level  and  pressure,  air 
ingestion  was  observed  to  begin  when  the  tank  pressure  at  the 
puncture  (vapor  space  pressure  plus  hydrostatic)  was  within 
0.25  KPa  (1  in.  of  water)  or  less  of  atmospheric  pressure; 
that  is: 

PT  '  patm  +  9»l  (ZL  -ZLh>  <  0-°°25  Von  (II1'2) 

Since  the  vapor  space  pressure  P^  could  not  be  held  exactly 
constant  during  a  test,  with  pressure  variations  of  +  1  cm  wa¬ 
ter  not  being  uncommon,  the  value  of  the  left  hand  side  of 
equation  (III. 2)  at  the  instant  air  ingestion  started  could 
not  be  calculated  with  exact  certainty,  but  it  was  in  the  range 

of  0  to  0.0025  times  P  .  . 

atm 


VACUUM 

RELIEF 

VALVE 


PRESSURE 

GAGE 


FIGURE  III. 2.  TEST  APPARATUS  FOR  AIR  AND  WATER 
INGESTION  STUDIES 


35 


The  next  set  of  tests  continued  the  venting  past  the 
point  of  air  ingestion  until  the  liquid  level  dropped  to  with¬ 
in  2  .  5  cm  (1  in.)  of  the  puncture  centerline.  These  tests 
confirmed  the  assumptions  used  in  developing  the  model  of  air- 
ingestion/intermittent  outflow;  that  is:  (1)  the  flow  does 
not  cease  when  "  patn  +  90^  ^ZL-2Lh^  -^rst  equals  zero,  if 
ZL  >  ZLh'  ^  air  ^-n9estec^  the  fozro  of  discrete  bubbles, 
comparable  in  average  diameter  to  about  two  or  three  times  the 
puncture  diameter;  (3)  during  the  time  the  bubble  blocks  the 
puncture,  the  outflow  ceases;  and  (4)  the  vapor  space  pressure 
increases  each  time  a  bubble  is  ingested.  Table  III. 2  sum¬ 
marizes  the  results  of  these  tests.  Again,  each  test  was  con¬ 
ducted  at  least  three  times,  with  excellent  repeatability  of 
the  results  (less  than  1  second  variation  in  discharge  time) . 

Rise  times  of  the  bubbles  were  measured  in  several  tests. 
The  bubbles  appeared  to  attain  their  terminal  velocity  very 
quickly;  the  average  velocity  was  computed  to  be  50.2  cm/sec 
+4.5  cm/sec. 

Water  Ingestion  Studies.  For  these  tests,  the  tank  was 
installed  in  a  shallow  but  larger  diameter  sheet-metal  tank, 
which  was  then  filled  with  water  until  the  puncture  was  sub¬ 
merged  to  a  preselected  depth.  Three  different  liquids  were 
tested:  carbon  tetrachloride,  which  is  heavier  than  water; 

isopropyl  alcohol,  which  is  lighter  than  water;  and  aqueous 
glycerine,  which  is  heavier  than  water.  Carbon  tetrachloride 
and  isopropyl  alcohol  are  both  virtually  immiscible  with  water, 
but  water  dissolves  readily  in  aqueous  glycerine.  For  these 
tests,  the  simulated  safety  relief  valve  was  closed  completely 
to_emphasize  the  water  ingestion  aspects  of  the  discharges. 

The  vapor  space  pressure  during  these  tests  dropped  almost 
immediately  to  a  vacuum  about  numerically  equal  to  the  hydro¬ 
static  head  of  the  liquid,  and  water  was  ingested  into  the 
tank  almost  immediately.  After  passing  through  the  puncture, 
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the  water  stream  broke  up  into  droplets  which  rose  to  the  sur¬ 
face  for  the  carbon  tetrachloride  tests,  or  sank  to  the  tank 
bottom  for  the  isopropyl  alcohol  tests,  or  mixed  with  the  tank 
liquid  for  the  aqueous  glycerine  tests.  The  water  ingestion 
continued  thereafter  until,  for  the  carbon  tetrachloride  and 
isopropyl  alcohol  tests,  the  level  of  water  in  the  tank  cov¬ 
ered  the  puncture.  For  the  aqueous  glycerine  tests,  the  den¬ 
sity  of  the  tank  liquid  continuously  decreased  as  a  result  of 
mixing  with  the  ingested  water,  and  the  outflow  decreased 
steadily;  however,  even  after  45  minutes,  some  outflow  was 
still  observed.  (The  aqueous  glycerine  tests  were  conducted 
to  study  the  basic  phenomena  of  ingestion  more  completely; 
however,  the  venting  rate  model  does  not  include  any  provisions 
for  mixing  of  the  ingested  water  with  the  cargo.)  Table  III. 3 
summarizes  the  results  of  the  tests  using  carbon  tetrachloride 
and  isopropyl  alcohol.  For  all  these  tests,  the  tank  was  sub¬ 
merged  in  water  to  a  depth  7.6  cm  (3  in.)  above  the  centerline 
of  the  puncture,  and  only  the  smaller  puncture  (2.2  cm)  was 
used. 

The  tank  vacuum  pressure  during  the  carbon  tetrachloride 
test  slowly  decreased  as  water  was  ingested,  with  the  net 
effect  being  that  the  total  pressure  (vapor  space  +  liquid 
head)  at  the  puncture  remained  fairly  constant.  At  the  in¬ 
stant  that  water  was  first  ingested,  for  example,  the  total 
pressure  at  the  puncture  was  11.4  cm  of  water,  or  just  slightly 
greater  than  the  hydrostatic  head  of  water  outside  the  tank 
(7.6  cm  of  water);  as  the  carbon  tetrachloride  reached  the  top 
cf  the  puncture,  a  point  where  the  water  level  in  the  tank  was 
27.6  cm  above  the  puncture,  the  total  pressure  at  the  puncture 
was  8.5  cm  of  water.  It  can  be  concluded  that  the  mass  flow 
rates  of  the  water  and  the  cargo  adjusted  themselves  to  main¬ 
tain  this  constancy  of  pressure,  and  the  analytical  model  of 
water  ingestion  for  a  cargo  heavier  than  water  is  based  upon 
this  conclusion. 


TABLE  III. 3.  WATER  INGESTION  STUDIES 


III. 3a  Carbon  Tetrachloride,  p  =  1.59  g/cm 


Tank  Pressure 

Time , 

Liquid  Height,  cm 

cm  of  H20  Vacuum 

sec 

27.9 

33.0 

mgm 

22.9 

30.5 

17.8 

27.9 

245 

12.7 

24.9 

380 

7.6 

22.1 

510 

2.5 

20.3 

635 

top  of  puncture 

19.  1 

665 

At  end  of  test,  water  height  in  tank  =  27.6  cm 


III. 3b  Isopropyl  Alcohol,  p  =  0.783  g/cm' 


Tank  Pressure 

Liquid  Height,  cm  cm  of  ^0  Vacuum 


Time , 
sec 


At  end  of  test,  water  height  in  tank  *  5.6  cm,  or 
level  with  puncture  top. 
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For  the  isopropyl  alcohol  tests,  the  vapor  space  pressure 
remained  nearly  constant,  as  did  the  level  of  the  isopropyl 
alcohol  in  the  tank.  The  ingested  water,  which  sank  to  the 
bottom  of  the  tank,  appeared  to  displace  the  tank  liquid  such 
that  the  inflow  and  outflow  volumetric  flow  rates  always  re¬ 
mained  equal,  and  neither  the  tank  pressure  nor  the  cargo 
liquid  height  varied  greatly.  In  effect,  all  the  cargo  below 
the  top  of  the  puncture  was  discharged  and  replaced  by  wauar. 
(The  slight  drop  of  0.6  cm  in  the  level  of  the  isopropyl  alco¬ 
hol  may  be  due  to  its  solubility  in  water,  which  is  small  but 
nonzero.)  The  analytical  model  for  water  ingestion  into  a 
tank  containing  a  cargo  lighter  than  water  is  based  upon  this 
conclusion  about  the  flow  rates. 

III. 2. 3  Test  Series  No.  3:  Discharge  of  Volatile 

Cargos 

The  final  series  of  tests  were  designed  to  evaluate  such 
thermodynamic  phenomena  for  volatile  cargos  as  liquid  strati¬ 
fication,  finite  evaporation  rates,  and  two-phase  flow,  and  to 
obtain  quantitative  data  for  verification  of  the  venting  rate 
model. 

A  closed  rectangular  tank  30.5  cm  by  30.5  cm  (12  in.  by 
12  in.)  and  61.0  cm  (24  in.)  high,  made  of  0.63  cm  (0.25  in.) 
thick  aluminum  plate,  with  a  total  capacity  of  56.6  liters 
(15  gallons),  was  used  for  these  tests;  see  Figure  III. 3.  The 
tank  exterior  was  insulated  by  fiberglass  5  cm  (2  in.)  thick 
to  approximate  adiabatic  test  conditions.  Two  600-watt  heat¬ 
ing  tapes  were  attached  to  the  wall  to  bring  the  simulated 
cargo  to  the  desired  temperature  before  a  test.  A  small  ex¬ 
terior  pump  continuously  circulated  the  liquid  in  the  tank 
during  the  pretest  period  to  prevent  any  initial  stratifica¬ 
tion;  the  spatial  variation  of  liquid  temperature  was  always 
kept  below  0.2'7oC  (0.5°F)  prior  to  a  test  by  this  procedure. 

Two  simulated  punctures,  0.32  cm  (0.125  in.)  and  2.2  cm 
(0.375  in.)  diameter,  were  drilled  into  one  wall  near  the  ver¬ 
tical  centerline  and  7.6  cm  (3  in.1  above  the  bottom;  most  of 
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FIGURE  III. 3.  INSULATED  AND  HEATED  TEST  TANK  FOR 
DISCHARGE  TESTS  OF  VOLATILE  LIQUIDS 
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the  testing  was  done  with  the  larger  puncture.  A  valve  mounted 
on  the  tank  top  was  used  to  vent  vapor  during  the  preliminary 
heating,  and  for  some  tests,  as  a  simulated  relief  valve. 

Tank  liquid  levels  were  measured  visually  from  a  graduated 
sight  tube,  marked  along  each  2.54  cm  (1.0  in.)  of  its  length; 
in  some  tests,  these  measurements  were  verified,  with  excel¬ 
lent  correlation,  by  thermocouple  measurements  as  explained 
below. 

The  tank  was  instrumented  with  20  Type  J  bayonet  thermo¬ 
couples  inserted  in  one  wall  through  drilled  holes  and  extend¬ 
ing  to  the  vertical  centerline  of  the  tank;  the  vertical  spac¬ 
ing  of  the  holes  was  2.54  cm.  The  vapor  space  pressure  was 
measured  with  a  diaphragm  type  of  transducer  (Validyne  DP7) . 

The  thermocouple  and  pressure  measurements  were  continuously 
recorded  on  a  light  beam,  strip  chart  oscillograph  (Honeywell 
Model  1858  Visicorder) .  A  stop  watch  was  used  to  measure  dis¬ 
charge  times;  these  measurements  were  verified  by  the  timing 
marks  on  the  strip  chart. 

Isopentane  and  methylene  chloride  (also  called  dichlcro- 
methane) ,  both  of  which  have  boiling  points  slightly  above 
room  temperature  at  atmospheric  pressure,  were  used  in  these 
tests  as  simulated  volatile  cargos.  The  properties  of  these 
liquids  are  shown  in  Table  III. 4.  (Thermodynamic  data  on  methy- 
methvlene  chloride  are  incomplete  and  somewhat  contradictory; 
thus,  these  tabulations  are  not  recommended  for  the  Chemical 
Properties  File  without  further  verification.) 

Test  Preparations.  Preliminary  testing  revealed  that  ob¬ 
taining  reproducible  data  for  venting  of  volatile  liquids 
would  be  difficult,  primarily  because  of  the  difficulty  in 
controlling  the  initial  air  content  of  the  vapor  space.  Any 
small  amount  of  air  mixed  with  the  vapor  in  the  vapor  space 
would  result  in  an  initial  tank  pressure  that  was  larger  than 
would  be  expected  on  the  basis  of  saturation  conditions  for  a 
pure  vapor.  This,  in  turn,  would  cause  a  rapid  decrease  in  the 
tank  pressure  during  venting.  That  is,  the  partial  pressure  of 


; 
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TABLE  III. 4.  SATURATED  THERMODYNAMIC  PROPERTIES 
OF  ISOPENTANE  AND  METHYLENE  CHLORIDE 

P  «  saturation  pressure  (KPa) ;  COT  ,  C  =  specific  heats  of  saturated 

V  rL  i  V 

liquid  and  vapor  (cal/g-°C);  X  =  latent  heat  of  evaporation  (cal/g); 

3 

p^  *  liquid  density  (g/cm  );  Z  =  vapor  compressibility;  and  T  *  tempera¬ 
ture  (°C)  . 

a.  ISOPENTANE  (molecular  weight  =»  72.147)  [14] 

log1Q  Pv  -  5.9666  -  1045 .87/ (T  +  236.18) 

PL  =  0.6405  -  1.00255  *  10~3(T) 

Z  *  0.984  -  1.126  x  10‘4(T)  -  2.182  *  10_6(X2) 

CpL  =  0.5218  +  1.0584  x  10“3(T) 

Cpv  =  0.360  +  2.12  x  10_4(T) 

2 

X  -  [65.2257  -  7.4638  x  10_2(T)  -  1.864  x  10_4(T2)]  (T  *  273  — - 

(T  +  236.13) ^ 

-  0.57530  Pv  (T  +  273.2)/pl  (T  +  236.18)2 

b.  METHYLENE  CHLORIDE  (molecular  weight  =  84.93)  [15] 

log10  Fy  =  5.4971  -  822.48/ (T  + 194.92) 

PL  =*  1.35  -  1.75  x  10'3(T) 

Z  -  0.9995  -  1.588  *  10_3(T)  +  8.505  x  10_6(T2) 

CpL  =  0.285  +  3.50  x  10_4(T) 

Cpv  -  0.129  -  1.651  x  10_3 (T) 

2 

-7  -4  2  (T  +  273  23 

X  =  [44.26  -  7.032  x  io  t  (T)  +  3.766  x  io  (T  )  ]  — — ~-~L— 

(T  +  194.92) 

-  0.45242  P  (T  +  273 .2) /pT  (T  +  194.92)2 

V  u 
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the  constant  mass  of  air  decreased  in  accordance  with  the  in¬ 
crease  in  vapor  space  volume  as  the  liquid  discharged.  The  par¬ 
tial  pressure  of  the  vapor,  on  the  other  hand,  remained  more 
nearly  constant  since  its  mass  was  increased  by  evaporation  from 
the  liquid  surface  at  a  rate  that  tended  to  compensate  for  the 
increase  in  vapor  space  volume .  Depending  upon  the  quantity  of 
air  initially  contained  in  the  vapor  space,  the  drop  in  the  vapor 
space  pressure  would  lead  to  air  ingestion  through  the  puncture 
fairly  quickly  for  some  tests,  while  for  other  tests,  with  seem¬ 
ingly  identical  conditions,  air  ingestion  would  occur  only  near 
the  end  of  the  discharge,  if  at  all.  After  this  was  realized, 
great  care  was  taken  to  eliminate  all  the  air  from  the  vapor 
space,  by  the  following  procedure. 

The  tank  was  filled  with  liquid  from  the  main  reservoir 
until  liquid  flowed  out  the  valve  in  the  top  of  the  tank.  The 
valve  was  then  closed,  the  filling  halted,  and  the  liquid 
heated  until  a  small  positive  pressure  existed.  The  valve 
was  then  opened  to  vent  a  liquid/vapor  mixture,  all  the  while 
continuing  to  heat  the  liquid  to  maintain  a  positive  pressure 
in  the  tank  and  thus  prevent  air  from  entering  through  the 
valve.  The  boil-off  of  liquid  was  continued  until  the  desired 
liquid  level  in  the  tank  was  obtained.  The  valve  was  then 
closed.  The  heating  and  the  circulating  of  the  liquid  through¬ 
out  the  tank  interior  by  the  pump  were  continued,  however,  un¬ 
til  the  desired  temperature  and  pressure  conditions  were  ob¬ 
tained.  Liquid  (or  vapor)  discharge  was  then  started  by  man¬ 
ually  removing  a  rubber  cork  from  the  simulated  puncture. 

This  procedure  usually  led  to  reproducible  tests,  although 
some  variation  was  still  encountered  occasionally,  apparently 
because  of  air  that  became  dissolved  in  the  liquid  isopentane 
and  methylene  chloride  in  their  reservoirs.  Altogether,  nearly 
40  tests  were  needed  to  obtain  reliable  venting  data  for  each 
of  the  desired  sets  of  initial  conditions. 

Liquid  Stratification.  Several  of  the  tests  were  de¬ 
signed  to  evaluate  the  assumption  that  venting  of  the  Liquid 
would  cause  only  a  negligible  degree  of  stratification. 
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Figures  III. 4  and  III. 5  show  typical  results  of  the  thermo¬ 
couple  measurements  for  fairly  rapid  discharges  of  isopentane 
and  methylene  chloride  through  the  2. 2 -cm  puncture.  Some 
stratification  is  evident  in  the  liquid  near  the  liquid-vapor 
interface  as  a  consequence  of  the  heat  transfer  needed  to  pro¬ 
vide  the  latent  heat  of  vaporization,  and  there  is  a  1°-2°C 
temperature  drop  from  the  liquid  to  vapor.  However,  this 
small  amount  of  stratification  has  a  negligible  effect  on  the 
energy  content  of  the  cargo,  and  the  assumption  that  the  tem¬ 
perature  is  uniform  is  acceptable  for  the  purposes  of  formu¬ 
lating  the  energy  balances,  as  will  be  shown  later  by  the  test 
and  model  comparisons.  For  slower  discharge  rates  through  the 
0.32-cm  puncture,  stratification  was  even  less  noticeable. 

The  thermocouple  strip-chart  records,  which  clearly  dis¬ 
played  the  1°-2°C  temperature  decrease  at  the  liquid/vapor 
interface,  provided  a  convenient  check  on  the  sight  tube  mea¬ 
surements  of  liquid  level.  The  comparison  of  these  two  inde¬ 
pendent  determinations  was  excellent. 

Evaporation  Rate.  Figures  III. 6  and  III. 7  give  typical 
results  of  vapor  space  pressure  and  liquid  level  for  liquid 
discharges,  for  the  same  two  tests  given  in  Figures  III. 4  and 
III. 5.  The  measured  vapor  space  pressure  decreased  rapidly 
upon  opening  the  puncture,  and  thereafter  decreased  more 
slowly.  Also  shown  is  the  calculated  saturation  pressure  of 
the  liquid,  based  on  the  average  liquid  temperature  estimated 
from  the  thermocouple  records  and  on  the  liquid  surface  tem¬ 
perature;  it  can  be  seen  that  these  two  pressures  differ  very 
little,  a  further  confirmation  of  the  negligible  effect  of 
stratification.  However,  the  actual  vapor  space  pressure  is 
significantly  less  than  the  saturation  pressure,  as  it  should 
be  for  the  reasons  discussed  in  Section  II. 3.  A  drop  in  vapor 
space  pressure  was  observed  even  for  the  much  slower  venting 
through  the  0.32-cm  puncture,  but  the  magnitude  of  the  drop 
was  much  smaller,  as  shown  for  a  typical  case  in  Figure  III. 8. 


FIGURE  III. 7.  DISCHARGE  OF  METHYLENE  CHLORIDE;  PUNCTURE  DIAMETER  =  2.2  cm 

TEST  NO.  18. 
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The  increase  in  the  mass  of  vapor  contained  in  the  vapor 
space  for  these  tests  was  computed  by  using  the  measured  vapor 
space  pressure  and  temperature  records  in  conjunction  with  the 
measured  liquid  levels  (i.e.,  vapor  space  volumes).  The  evap¬ 
oration  rate  model  described  earlier  in  Section  II. 3  was  also 

used  to  predict  the  evaporated  mass,  since  W  could  be 

6  f  sat. 

computed  from  the  temperature  measurements  of  the  liquid, 
using  the  known  saturation  pressure  vs  temperature  relation 
for  the  tested  liquid  to  compute  the  needed  saturation  pres¬ 
sure.  Comparisons  of  the  model  predictions  (integrated  with 
respect  to  time  to  give  vapor  mass  change  rather  than  evapora¬ 
tion  rate)  to  the  tests  are  shown  in  Figure  III. 9.  As  can  be 
seen,  the  comparison  is  very  close.  These  kinds  of  test  and 
model  results  were  used,  in  fact,  to  deduce  the  best  value  of 
the  empirical  constant  (1.9)  in  equation  (11.12).  Predictions 
of  the  model  given  by  equation  (II. 11)  are  also  shown;  these 
comparisons  are  much  poorer. 

It  ought  to  be  mentioned  that  there  is  some  uncertainty 
in  the  accuracy  of  the  adiabatic  behavior  of  the  liquid  late 
in  each  test.  The  large  thermal  inertia  of  the  test  tank 
probably  moderated  the  decrease  of  the  liquid  temperature 
somewhat,  so  that  both  the  true  "adiabatic"  liquid  tempera¬ 
ture  and,  therefore,  the  saturation  pressure  should  actuallv 
be  slightly  less  near  the  end  of  the  tests  than  are  shown  in 
Figures  III. 4  through  III. 8. 

Two-Phase  Flow.  The  tank  liquid  level  data  as  a  function 
of  time  were  used  both  to  evaluate  the  implicit  assumption 
made  in  AMSHAH  [1]  that  two-phase  flow  effects  were  negligible, 
and  to  determine  the  best  correlating  value  of  the  empirical 
factor  8  in  the  present  two-phase  flow  model. 

For  the  test  results  shown  in  Figure  III. 8,  the  average 
discharge  rate  calculated  from  the  test  data  is  23.2  g/sec 
for  the  130  seconds  required  for  the  liquid  level  to  drop 
2.54  cm  from  an  average  level  of  34.3  cm  above  the  puncture 
centerline.  Since  the  tank  gage  pressure  at  this  time  was 
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TIME  AFTER  PUNCTURE  OPENED  (SEC) 


FIGURE  III. 9.  COMPARISON  OF  PREDICTED  AND  MEASURED  EVAPORATION 

ISOPENTANE  TEST  NO.  6. 


36.3  err.  of  water,  the  flow  model  used  in  [1]  predicts  that  the 
discharge  rate  should  be 


or  35.5  CD  g/sec.  The  discharge  coefficient  required  by  this 
model  to  acres  with  the  measured  W„  is  0.65,  but  the  actual 

'  uO 

(as  determined  by  prior  testing)  is  0.73.  The  model  there¬ 
fore  overpredicts  the  test  results  by  12%.  Likewise,  for  the 
test  data  shown  in  Figure  III. 4,  the  single-phase  flow  model 
predicts  a  discharge  rate  of  708  g/sec  (with  =  0.73)  at 
the  instant  when  the  tank  pressure  is  30.5  cm  F^O  and  the  liq¬ 
uid  level  is  35.5  cm,  although  the  discharge  rate  calculated 
from  the  test  data  is  about  600  g/sec;  thus,  the  model  also 
overpredicts  this  test  by  about  16%. 

All  the  tests  were  analyzed  at  representative  points  to 
determine  the  best  correlating  value  of  3  for  the  two-phase 
flow  model  (equations  (II. 7a)  and  (II.7b)).  The  best  average 
value  of  3  was  found  to  be  0.12,  and  there  was  little  test-to- 
test  variation.  For  example,  the  discharge  rates  predicted  by 
the  model  with  3  =  0.12  for  the  two  tests  described  above  are 
23.7  g/sec  and  636  g/sec  ,  respectively,  both  in  excellent 
agreement  with  the  tests. 

★  _  _ 

For  the  methylene  chloride  test  shown  m  Figure  II-.3, 

the  test  measurements  at  the  indicated  time  were  ?T  =  101.60 

KPa ,  P  =98.0  KPa,  T  *  41.5’C,  and  P„  -  =  -7.5  cm  H-,0  = 

a  tn  v 

-0.75  KPa.  The  corresponding  calculated  liquid  properties 
are  o.  =  1.278  g/cm3  ,  Tg  =  39.68°C,  =  74.6  cal/g  ,  v^s  = 

300.0  cm3/g  ,  and  C?T  (average)  =  0.299  cal/g-°C  .  From  eq. 
(II. 7b)  with  3  =  0.12,  vT^  =  1.036  cm3/g  (or  pT  =  0.965 
g/cm  )  .  The  discharge  velocity  (the  radical  in  eq.  (II. 7a) 
is  425  cm/sec.  Thus,  =  0.07917  *  0.73  *  425/1.036  =  23.7 

(Footnote  continued  on  next  page) 
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Vapor  Discharge.  Several  tests  were  conducted  in  which 
only  vapor  was  discharged,  using  the  valve  in  the  top  of  the 
tank  as  a  simulated  puncture.  There  was  an  initial  blowdown 
of  the  tank  pressure  to  very  near  atmospheric  pressure,  fol¬ 
lowed  by  a  much  slower  discharge  as  the  liquid  remaining  in 
the  tank  evaporated.  These  tests  were  not  analyzed  in  detail 
for  two  reasons:  test  durations  were  so  long  that  the  tank 
processes  could  not  validly  be  considered  as  adiabatic,  and 
the  test  pressure  instrumentation,  which  was  designed  primar¬ 
ily  for  liquid  discharges,  could  not  resolve  the  small  pres¬ 
sure  difference  between  the  tank  vapor  space  and  the  atmosphere 
which  would  be  needed  to  make  discharge  rate  comparisons  with 
the  model.  It  was  noted,  however,  that  the  temperature  of  the 
vapor  space  was  consistently  much  lower  than  the  liquid  as  a 
result  of  the  initial  blowdown;  a  typical  temperature  differ¬ 
ence  was  17 °C. 

Simultaneous  Liquid  and  Vapor  Discharge.  No  unusual  or 
'unexpected  phenomena  were  encountered  during  these  tests.  The 
vapor  space  pressure  temperature  again  was  nearly  17 °C  less 
than  the  liquid.  These  tests  were  not  analyzed  in  detail  for 
the  same  reasons  given  for  the  vapor  discharge  tests. 

Air  Ingestion.  In  a  few  tests,  air  was  trapped  in  the 
vapor  space  deliberately  during  the  tank  filling  process.  In 
this  way,  a  partial  vacuum  could  be  caused  almost  immediately 
after  the  start  of  the  liquid  discharge,  and  air  was  ingested 
through  the  puncture  during  most  of  the  test.  The  air  inges¬ 
tion  began  at  a  tank  pressure  and  liquid  level  that  was  in 
agreement  with  equation  (III. 2) .  No  unusual  or  unexpected  phe¬ 
nomena  were  observed,  so  the  tests  were  not  analyzed  further. 

For  the  isopentane  test  shown  in  Fiaure  III. 4,  =  100.73 

i. 

K?a ,  P=.  =  97.74  K?a,  T  =  2S.9°C,  ?„  -  P,„  =  -42.0  cm  H0C  = 

-4.11  KPa ,  c L  =  0.611  g/cm  ;  Tg  =  26.83°C,  Xg  =  84.7  cal/g  , 
v..  =  346  cm-Vg  ,  and  C_.  (averace)  =  0.551  cal/g -°C.  The 
discharge  specific  volume  =  2.181  cm  /g  and  the  discharge 
velocity  =  490  cm/s ec .  Thus,  W  =  3.879  *  0.73  *  490/2.131 

LO 

=  636  g/sec. 
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IV.  VALIDATION  OF  MODEL 

3y  using  results  of  typical  tests,  all  the  empirical  con¬ 
stants  in  the  venting  rate  model  have  been  determined  by  the 
methods  discussed  in  Sections  II  and  III.  The  predictive  re¬ 
liability  of  the  overall  model  and  its  computerized  solution 
are  demonstrated  in  this  section  by  comparison  to  results  of 
all  the  test  series. 

IV. 1  Discharge  of  Nonvolatile  Liquids  into  Air; 

Operable  Relief  Valve 

The  discharge  of  a  nonvolatile  liquid  into  air  was  stud¬ 
ied  in  Test  Series  No.  1  as  a  means  of  determining  the  varia¬ 
tion  of  discharge  coefficients  with  the  shape  of  the  puncture. 
With  supplied  as  input  to  the  computerized  model,  the 
total  time  required  to  discharge  the  cargo  initially  above  the 
puncture  centerline  was  computed.  The  predicted  venting  times 
were  virtually  identical  to  the  test  results.  For  example, 
the  predicted  time  was  66.5  seconds  to  discharge  enough  liquid 
to  drop  the  level  from  64.8  cm  to  14.6  cm  above  the  puncture 
centerline  for  the  tank  shown  previously  in  Figure  II. 1.  The 
average  discharge  time  for  a  set  of  three  tests  was  also  66.5 
seconds . 

IV. 2  Discharge  of  Nonvolatile  Liquids:  Air  Ingestion 

Experimental  results  for  liquid  discharges  in  which  air 
was  ingested  (the  vacuum  relief  valve  was  jammed)  were  pre¬ 
sented  in  Table  II..  .2.  All  these  data  were  analyzed  in  con¬ 
junction  with  predictions  of  the  venting  rate  model  to  deter¬ 
mine  the  best  correlating  value  of  the  ingested  air  bubble 
volume  factor.  This  value  turned  out  to  be  5.6;  the  implied 
average  bubble  diameter  is  therefore  about  1.3  times  the  aver¬ 
age  open  dimension  of  the  puncture,  which  is  in  good  agreement 
with  the  visual  observations  of  bubble  size. 


Test  Series  No. 


i  (wnere 


Comparisons  of  a  typical  test  of 
the  simulated  cargo  was  water)  and  the  predictions  of  the 
model  are  given  below.  In  the  test,  the  simulated  relief 
valve  was  manually  adjusted  to  maintain  a  pressure  in  the  va¬ 
por  space  vacuum  of  15.25  cm  of  water,  so  the  relief  set-point 
used  in  the  prediction  was  also  input  as  this  value  (i.e., 
1.494  K?a) . 


INITIAL  HEIGHT 
ABOVE  PUNCTURE, 
cm 

TIME  TC 

INGESTION  START , 
sec 

HEIGHT  AT 
INGESTION  START, 
cm 

TOTAL 

OISCHARGE  TIME, 
sec 

Test 

Modal 

Test 

Model 

B&9 

Model 

2C.7 

26.7 

21.6 

1 

22.3 

17.8 

21. 3 

88.1 

36.0 

As  can  be  seen,  the  model-tc-test  comparison  is  excel¬ 
lent,  both  with  respect  to  when  ingestion  begins  and  to  the 
total  time  required  to  discharge  the  liquid. 

IV. ^  Discharge  of  Nonvolatile  Liquids:  Water  Ingestion 

Experimental  results  for  the  discharge  of  liquids  when 
water  was  ingested  were  given  in  Table  III. 3.  These  tests 
were  analyzed  in  conjunction  with  the  venting  rate  model  to 
determine  the  best  correlating  value  of  the  empirical  parame¬ 
ter  K  of  the  model  (see,  for  example,  equations  (11.17));  the 
best  overall  correlation  for  all  the  tests  was  found  to  be 
given  for  K  =  83. 

When  the  simulated  cargo  was  heavier  than  water,  the  model 
predicted  the  tests  excellently,  as  shown  in  the  table  in  the 
next  page.  In  the  tests,  the  relief  valve  was  kept  closed  to 
create  a  vacuum  pressure  such  that  water  was  ingested  almost 
immediately  upon  starting  the  discharge.  In  the  model,  nearly 
the  same  vacuum  was  predicted,  with  water  being  ingested  after 
0.8  second.  For  both  the  test  and  the  model,  the  external  water 
level  was  12.7  cm  above  the  puncture  centerline.  As  can  be  seen, 
the  model  and  the  test  are  well  correlated,  including  the  trend 
and  magnitude  of  the  vapor  space  vacuum  as  a  function  of  time. 


DO 


[■■■SlfeHH 

cm  H,0 

FINAL  VACUUM, 
cm  K.,0 

TOTAL 

DISCHARGE  TIME, 
sec 

Test 

Model 

Tes- 

Model 

15HB 

■1 

Test 

Model 

27.9 

27.9 

33.0 

28.9 

17.8 

13.1 

655 

641 

3y  referring  to  equations  (11.17),  it  can  be  seen  that  the 

model  prediction  of  the  discharge  rate  depends  on  P  .  and  Z  . , 

XX  Li 

the  values  of  tank  pressure  and  liquid  height  at  the  time  in¬ 
crement  just  before  water  ingestion  starts.  For  a  given  set 
of  initial  conditions,  the  numerical  incrementing  scheme  will 
always  arrive  at  these  same  values  of  P  .  and  Z  . ,  and  will 

T X  LX 

oredict  the  same  values  of  P„  _  and  ZT  „  .  For  sliantlv 

i ,Eq.  L,iC. 

different  initial  conditions,  ZT .  and  P_.  would  also  be  slightly 

LtX  rx 

different.  On  a  percentage  basis,  however,  the  changes  in 
-  P^  __  and  Z. •  -  ZT  _  would  be  substantial.  Thus,  the 
discharge  rate  prediction  might  change  significantly  for,  say, 
a  small  change  in  the  initial  liquid  level,  although  the  ac¬ 
tual  discharge  rate  would  be  expected  to  change  only  slightly. 
Consequently,  X  would  also  need  to  be  changed  to  bring  the 
predictions  into  agreement,  which  is  an  undesirable  result. 

To  determine  whether  the  predictions  are  in  fact  sensitive  to 
small  changes  in  the  initial  conditions,  the  initial  liquid 
level  in  the  computations  was  decreased  by  0.1  cm,  to  27.3  cm 
above  the  puncture.  The  corresponding  change  in  ZT .  -  Z_  _ 

^jX  la  f  X C7  . 

was  fifty  percent.  However,  the  predicted  total  discharge  time 
changed  by  only  3  seconds,  to  538  seconds,  and  the  changes  in 
the  computed  initial  and  final  vacuum  pressures  were  also  neg¬ 
ligible.  The  reason  for  the  lack  of  sensitivity  is  that  the 
slight  increase  of  the  initial  vapor  space  volume  resulted  in 
a  slight  increase  of  Pm.  ,  which  increased  P„.  -  P^,  _  just 
enough  to  cancel  cut  the  fifty  percent  decrease  of  ZT .  -  Z,  „ 

la  X  1j  /  H  q  • 

It  is  concluded  that  there  is  very  little  sensitivity  to  slight 
changes  in  the  initial  conditions.  The  corollary  to  this  con¬ 
clusion  is  that  the  value  X  =  S3  is  general  and  not  specific 
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to  the  particular  values  of  ZT^  and  Pm.  that  result  from  a 
given  set  of  initial  conditions  or  the  numerical  solution  tech¬ 
nique  . 

Correlation  of  the  model  to  the  test  results  for  cargos 
lighter  than  water  was  not  as  good.  For  example,  for  the 
isopropyl  alcohol  tests  described  in  Section  III,  the  mea¬ 
sured  time  required  for  the  ingested  water  level  to  reach  the 
top  of  the  puncture  was  613  seconds,  at  which  time  the  dis¬ 
charge  ceased.  The  computed  results  were  that  194  seconds  are 
required  to  ingest  enough  water  to  cover  the  puncture;  more¬ 
over,  this  prediction  was  not  very  sensitive  at  all  to  the 
value  of  K  used  in  the  range  of  83  to  1000  ,  so  the  value  K  =  33 
selected  from  the  tests  with  cargos  heavier  than  water  is  also 
recommended  for  cargos  lighter  than  water.  The  primary  reason 
for  the  discrepancy  in  the  model  predictions  is  believed  to  be 
as  follows:  the  conditions  for  water  ingestion  were  met  at 
the  end  of  the  first  mass  discharge  incrementation  cycle  of 
the  model.  This  meant  that  PTi  (the  pressure  at  the  iteration 
preceding  the  start  of  ingestion)  was  selected  as  the  initial 
vapor  space  pressure;  that  is,  PT^  was  equal  to  atmospheric 
pressure.  This  value  of  tank  pressure  is  the  pressure  that  is 
subsequently  maintained  constant  by  the  water  ingestion  model , 
rather  than  the  12.7  cm  of  water  vacuum  observed  in  the  tests. 
The  result  is  that  the  predicted  discharge  is  too  large.  For 
a  ship  tank,  the  static  head  of  the  cargo  above  the  puncture 
would  be  much  larger  than  was  used  in  these  tests.  Thus,  the 
model  would  predict  that  a  definite  volume  of  liquid  would  be 
discharged  before  a  vacuum  large  enough  to  cause  water  inges¬ 
tion  would  be  created.  Therefore,  the  difficulty  encountered 
in  these  predictions  of  the  discharge  time  for  the  small  tanks 
used  in  the  present  tests  is  probably  not  a  serious  one. 

IV. 4  Discharge  of  Volatile  Liquids 

Experimental  results  for  the  discharge  of  a  volatile  liq¬ 
uid  were  discussed  in  Section  III.  3y  analyzing  the  test 


results,  an  evaporation  rate  model,  as  described  in  Section 
II. 3,  was  formulated  and  the  empirical  constants  of  the  model 
were  determined.  The  empirical  factor,  3,  in  the  model  of  two- 
phase  flow  was  determined  similarly.  The  entire  venting  rate 
model  for  volatile  cargos  is  now  compared  to  the  tests. 

Figure  IV. 1  shows  a  comparison  between  the  predictions  of 
liquid  height,  evaporated  mass,  and  vapor  space  pressure,  and 
representative  results  of  a  typical  isopentane  test  taken  from 
the  strip-chart  recordings.  The  predicted  total  time  to  dis¬ 
charge  the  cargo  is  about  63  seconds,  compared  to  about  67  in 
the  test;  equally  good  agreement  is  apparent  for  the  other  pa¬ 
rameters  of  interest.  In  particular,  the  vapor  space  pressure 
decreased  with  time  in  the  same  way  as  did  the  test;  this  indi¬ 
cates  that  the  evaporation  rate  model  includes  the  appropriate 
physical  phenomena,  even  though  the  total  mass  evaporated  into 
the  vapor  space  is  slightly  overpredicted.  Note  that  PT  ~ ?atm 
*  0  over  most  of  the  test  and  that  psat  always  remains  notice¬ 
ably  larger  than  patm*  Consequently,  basing  the  predictions 
on  ?sat  (as  AMSHAH  does)  would  overpredict  the  discharge  rate, 
especially  when  the  hydrostatic  head  was  small,  and  might  cause 
the  onset  of  air  or  water  ingestion,  when  they  occur,  to  be 
missed  entirely. 

The  discharge  rate  near  the  start  of  the  test  is  slightly 
overpredicted  because  of  the  slight  overprediction  of  the  evap- 
evaporation  rate,  and  this  causes  the  prediction  of  tank  liquid 
height  to  be  less  than  the  measured  height  throughout  the  rest 
of  the  discharge  time.  Note,  however,  that  the  discharge  rate 
(i.e.,  the  slope  of  the  liquid  height  curve)  is  well  predicted 
except  for  this  initial  discrepancy.  Air  ingestion  was  pre¬ 
dicted  to  occur  during  the  last  12  seconds  of  the  discharge. 

The  time  at  which  air  was  first  ingested  in  the  test  was  dif¬ 
ficult  to  determine  since  the  tank  was  not  transparent,  but  os¬ 
cillations  in  the  sight  tube  liquid  level  and  in  the  pressure 
record  became  noticeable  during  the  last  6  to  3  seconds,  which 
indicate  that  air  was  being  ingested  for  at  least  this  period. 
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FIGURE  IV. 1.  COMPARISON  OF  VENTING  RATE  MODEL  AND  TEST  RESULTS  FOR 
VENTING  OF  LIQUID  ISOPENTANE;  TEST  NO.  6 
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Thus,  the  model  and  test  are  also  in  good  agreement  concerning 
air  ingestion. 

Figure  IV. 2  shows  similar  comparisons  for  a  methylene 
chloride  test.  The  predicted  discharge  time  of  93  seconds 
compares  well  to  the  measured  result  of  nearly  100  seconds. 
Predictions  of  the  other  phenomena  also  compare  well,  except 
that  the  magnitude  of  the  initial  rapid  decrease  of  the  vapor 
space  pressure  was  underpradicted.  This  is  perhaps  the  re¬ 
sult  in  the  test  of  the  existence  of  a  small  quantity  of  air 
in  the  vapor  space . 

The  initial  temperature  of  the  liquid  for  the  predictions 
shown  in  both  Figures  IV. 1  and  IV. 2  had  to  be  adjusted  slightly 
(less  than  0.5aC)  in  comparison  to  the  test  measurement  in  or¬ 
der  to  bring  the  initial  vapor  space  pressure  into  agreement 
with  the  tests.  The  experimental  measurements  of  temperature, 
vapor  space  gage  pressure,  and  local  atmospheric  pressure  (all 
of  which  enter  into  the  comparison  of  predicted  and  measured 
gage  pressures)  could  have  been  in  error  by  enough  to  explain 
this  adjustment  of  initial  temperature,  but  it  is  also  pos¬ 
sible  that  the  industrial  quality  liquids  used  in  the  tests 
may  not  have  had  exactly  the  same  vapor  pressure  versus  tem¬ 
perature  relation  that  was  used  in  the  model  predictions. 

IV. 5  Other  Validation  Checks 

Data  from  other  tests,  such  as  venting  of  vapor  of  vola¬ 
tile  cargoes,  were  not  reliable  enough  to  serve  as  checks  on 
the  validity  of  the  model.  The  computer  model  has  been  exer¬ 
cised,  however,  for  all  the  various  kinds  of  possible  dis¬ 
charges  of  liquid  and  gas.  These  results  have  been  verified 
against  the  analytical  model  by  (1)  single-point  hand  calcula¬ 
tions  and  (2)  reasonableness  of  predictions.  As  far  as  could 
be  determined  by  these  limited  checks,  the  computerized  model 
accurately  reproduces  these  kinds  of  discharges. 
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V.  CONCLUSIONS 


The  Venting  Rate  Model  previously  developed  by  others  for 
the  Hazards  Assessment  Computer  System  has  been  completely  re¬ 
formulated.  All  errors  in  the  basic  equations  have  been  cor¬ 
rected,  and  the  limitations  imposed  by  inappropriate  assump¬ 
tions  have  been  removed.  The  revised  model  now  covers  nearly 
every  practical  combination  of  puncture  location,  cargo  proper¬ 
ties,  intermittent  outflows,  and  relief  valve  actuations.  In 
summary,  the  revised  model  incorporates  the  following  new  fea¬ 
tures  : 


o  punctures  above  or  below  the  waterline; 

o  relief  valve  actuations; 

o  the  possibility  of  air  in  the  vapor  space; 
o  simultaneous  liquid  and  gas  discharges; 
o  phase  changes  during  the  discharge  of  volatile 
cargos ; 

o  choked  flow  of  gases  and  volatile  liquids; 
o  air  or  water  ingestion;  and 

o  realistic  effects  of  evaporation  within  the  tank. 


Four  of  the  submodels  of  the  Venting  Rate  Model  contain 
empirical  constants  whose  "best"  values  have  been  determined  by 
comparison  to  experimental  data.  Three  of  the  constants  are 
related  to  observable  physical  phenomena — the  air  bubble  size 
factor  of  5.6  in  the  air  ingestion  model,  the  constant  1.9  in 
the  tank  cargo  evaporation  model,  and  the  constant  5  =  3.12  in 
the  two-phase  flow  model.  The  fourth  constant,  X  =  33,  is  used 
in  the  computerized  solution  of  the  water  ingestion  model.  The 
relative  importance  of  these  constants  in  the  predictions  de¬ 
pends  upon  the  type  of  cargo  discharge  encountered.  In  some 
cases  the  predictions  are  completely  independent  of  the  values 
of  the  constants.  For  example,  when  the  relief  valve  is  oper¬ 
able,  air  is  admitted  as  needed  to  keep  the  vapor  space  pres¬ 
sure  near  atmospheric.  Thus,  air  will  not  be  ingested  through 


the  puncture  until  the  liquid  discharge  is  nearly  finished, 
and  as  a  result  the  value  of  the  bubble  size  factor  is  unim¬ 
portant  in  the  model  predictions.  Further,  the  vapor  space 
pressure  is  controlled  by  the  relief  valve  and  not  by  evapora¬ 
tion  into  the  vapor  space,  so  the  numerical  constant  in  the 
evaporation  model  is  likewise  practically  irrelevant.  On  the 
other  hand,  water  can  still  be  ingested,  but  only  if  the  punc¬ 
ture  is  far  enough  below  the  waterline  to  overcome  the  hydro¬ 
static  head  of  the  cargo.  The  onset  of  water  ingestion  will 
be  considerably  delayed,  however,  compared  to  a  case  when  the 
relief  valve  is  jammed  shut.  The  effects  of  two-phase  flow, 
and  thus  the  value  of  S ,  are  only  important  for  discharges  of 
cargos  whose  vapor  pressure  is  greater  than  atmospheric  and 
are  carried  in  pressurized  tanks.  Otherwise  the  liquid  dis¬ 
charge  is  predicted  to  be  nearly  incompressible,  in  agreement 
with  observations.  In  any  event,  the  two-phase  flow  model  and 
the  tank  cargo  evaporation  model  affect  only  the  predicted  time 
required  for  the  discharge  and  not  the  total  quantity  of  liquid 
cargo  discharged.  A  faster  discharge  rate  and  a  more  conserva¬ 
tive  estimate  of  the  hazard  can  be  obtained  by  setting  either 
or  both  of  the  empirical  constants  (3  and  1.9)  in  these  models 
to  zero.  The  air  and  water  ingestion  models  do  influence  the 
computation  of  the  total  liquid  discharged,  but  the  value  of 
their  empirical  constants  (the  bubble  size  factor  and  K  =  83) 
only  change  the  discharge  rate.  A  conservative  estimate  of 
the  hazard  can  be  obtained  by  decreasing  the  value  of  K  (now 
equal  to  83)  in  the  water  ingestion  model  or  by  increasing  the 
value  of  the  bubble  size  factor  (now  equal  to  5.6)  in  the  air 
ingestion  model;  either  of  these  changes  will  increase  the  pre¬ 
dicted  discharge  rate. 

Predictions  of  the  computerized  Venting  Rate  Model  have 
been  compared  to  scale-model  test  results  for  many  kinds  of 
discharges.  Generally  good  correlation  of  the  model  to  the 
tests  was  found.  In  one  sense,  these  good  correlations  do  not 
necessarily  represent  a  true  verification  of  the  model  because 


part  of  the  data  used  in  making  the  comparisons  was  used  to 
determine  the  empirical  constants  of  the  various  submodels. 
Usually,  model  verification  requires  a  comparison  of  the  model 
to  an  independent  set  of  data  taken  under  different  experimen¬ 
tal  conditions  from  those  which  are  used  to  determine  the  em¬ 
pirical  constants.  3ecause  of  the  constraints  in  the  scope  of 
the  project,  it  was  not  possible  to  conduct  a  separate  set  of 
experiments  for  the  purpose  of  model  verification.  Therefore, 
an  alternative  technique  was  used. 

Each  submodel  containing  an  empirical  constant  was  first 
fitted  to  the  relevant  test  data  to'determine  the  constant. 

For  example,  in  determining  3  of  the  two-phase  flow  model,  test 
measurements  of  liquid  temperature,  liquid  level,  vapor  space 
pressure,  etc.,  were  used  to  make  predictions  of  mass  discharge 
rates,  using  the  flow  model  alone;  these  were  then  compared  to 
the  measured  discharge  rates  alone.  Note  that  actual  test 
data  were  used  to  determine  S  rather  than  predictions  of  liq¬ 
uid  temperature,  vapor  space  pressure,  and  the  other  required 
submodel  inputs  made  from  exercising  the  entire  Venting  Fate 
Model.  Since  only  the  flow  submodel  was  used  to  find  3,  it 
was  considered  an  appropriate  test  of  the  entire  Venting  Rate 
Model  to  input  the  initial  conditions  of  the  same  experiments 
to  compare  its  output  to  the  actual  test  results.  The  good 
correlation  of  the  entire  model  to  test  data  thus  does  show 
that  the  energy  balances  have  been  formulated  correctly,  that 
the  interaction  of  the  energy  balance  with  tank  cargo  evapora¬ 
tion  has  been  modeled  realistically,  that  the  interaction  of 
the  air  ingestion  model  with  the  vapor  space  pressure  model  has 
also  been  incorporated  realistically,  and  so  on.  It  can  be  con¬ 
cluded,  then,  that  the  Venting  Rate  Model  has  been  verified  to 
a  significant  degree.  However,  we  would  recommend  that  further 
experimentation  and  subsequent  model  validation  be  carried  out. 

The  Chemical  Properties  File  has  been  reviewed  and  found 
to  be  satisfactory  for  nonvolatile  cargos.  For  volatile  car¬ 
gos,  however,  several  additional  correlations  are  required: 


the  compressibility  and  specific  heat  of  saturated  vapor,  and 
the  latent  heat  of  evaporation.  It  also  appears  that  the  accu¬ 
racy  of  the  saturation  pressure  correlations  should  be  improved. 
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APPENDIX  A.  DERIVATION  OF  EQUATIONS 


A. 1  Tank  Energy  3alance 

The  tank  contents  are  assumed  to  be  "well  stirred,"  and 
the  temperature  of  the  gas  and  the  liquid  are  assumed  to  be 
equal.  Considering  a  control  volume  around  the  liquid  shown 
in  Figure  II. 1,  the  general  well-stirred  reactor  energy  equa¬ 
tion  is  : 
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The  work  term  W  includes  the  P^d V  contribution  of  the  moving 
boundaries.  Equation  (A.l)  is  specialized  for  adiabatic  pro¬ 
cesses,  Qcv  =  0.  The  internal  energy  of  the  control  volume  is: 


Ecv  -  h.  ("L  +  T  Vl  *  *=0  <A-2' 

where  u  is  the  liquid  internal  energy,  V  is  the  average  liq- 

Li  L» 

uid  velocity  in  the  tank  and  Zr  is  the  vertical  location  of 

“cm 

the  liquid  center  of  mass. 

Stagnation  enthalpies  are  defined  as  follows: 
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VGo 

(A. 3b) 

ha  = 

h,  +  \  Vr 

Ao  2 

(A. 3c) 

H 

TO 

hA + 1  vi 

(A. 3d) 
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where  the  first  two  relations  follow  from  the  fact  that  the 
liquid  (or  the  vapor)  in  the  tank,  and  in  the  discharge  stream 
both  have  the  same  reservoir  conditions.  Since  u  =  h. 
equation  (A.l)  can  be  re-written  as: 


±  (-\Hl  +  gZL_)  -  £r  (P TVr)  =  -Amp 


dz. 


^cm  dt  T '  L 


T*  T  dt 


We  (Bv  *  I  Ve  -  r  VG  *  5  st)  -  wLo  (H.  *  5  (A. 4) 


An  energy  equation  for  a  control  volume  around  the  vapor 
space  can  be  written  similarly: 


as  1  Vv  *  maha  *  ima  *  V  9  ZGcJ  -  gr  <  VV  - 


dZ, 


T  /  12  12 

Vt  +  We  (HV  +  I  ve  7  VG  +  ?  ZL 


-  *Vo  (HV  +  *  ZGh>  '  WAo  fHA  +  5  W 


)  +  WA  (h.  +  g  z  j 

(A. 5) 


Here,  the  assumption  has  been  made  that  the  sum  of  the  air  and 
the  vapor  partial  pressures  is  equal  to  the  tank  pressure,  so 
that  ?AhAvA  +  PvMvvv  can  be  set  equal  to  PAVG  +  pvvq  = 

(PA  +  V  VG  ”  PTVG* 

Equations  (A. 4)  and  (A. 5)  can  be  combined  by  using  conser¬ 
vation  of  mass  (e.g.,  equations  (II. 2))  to  eliminate  the  dM^/dt 
and  dM  /dt  terms,  and  by  using  the  definition  of  the  latent 
heat  of  vaporization: 

Kv  =  Hl  -  X  +  (V*  -  v\)/2  (A. 6) 


The  result  is: 
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dH 


L 


dH. 


V 


dH, 


dt 


"V  dt 


+  M 


A  dt 


dPT 

dF  +  WA  (H»  -  “a1 


*  +  I  Ve  +  (VG  " 


d.^ 


+  g  <  W.Z_  -  W_  ZT,  -  (W,. 

’  >  AT  Lo  Lh  Vo 


-  Zt  l\ZLCK  *  (ma  *  -V 


dt"  +  ^Lo 


WAoJ  ZGh 


"Gen 


1 


(A. 7) 


According  to  the  numerical  examples  given  in  [2] ,  both  the  ki¬ 
netic  and  potential  energy  terms  are  negligible  during  the  dis¬ 
charge  of  a  volatile  substance.  Thus,  since  dH/dt  =  dh/dt  = 

Cp  dT/dt  when  kinetic  energy  is  negligible,  equation  (A. 7)  can 
be  finally  written  as 


(MLCPL  +  MVCPV  +  MACPA)  dt 


WACPA 


(Tair  *  T) 


(A. 8) 


which  is  equation  (II.  1)  of  the  report. 

A. 2  Flow  Equations  for  Volatile  Substances 

Considering  an  isentropic  discharge  of  liquid  from  a  sat¬ 
urated,  or  equilibrium,  initial  state,  to  an  exit  pressure  P,, 
corresponding  to  a  saturation  temperature  Tg ,  where  the  fluid 
is  a  liquid-vapor  mixture,  the  thermodynamic  quality  of  the 
exit  mixture  is: 


SQUt  Sk=°  _  SL  S^°° 

SVS  -  sLs  *  sVs  -  sLg 


(A. 9) 


where  s.  is  the  entropy  of  the  saturated  liquid  in  the  tank. 
Li 

(s_  .  )  =  s.  =  sT  for  an  isentrooic  Drocess)  .  But 
out  in  L  *  - 


Ls 


/ 

Toa 


C?LdT 

T 


3. 


(A. 10a) 


2 


and 


s 


J  s 


3Ls 


i  /T 

3  S 


(A.ICb 


Thus,  the  enthalpy  change  is: 


in 


out 


=  h.  -  [h„  +  xA  ] 

in  Ls  s 


(A. 11) 


If  the  liquid  is  compressed,  i.e.,  if  the  tank  pressure  is 
greater  than  the  saturation  pressure,  there  is  an  additional 
enthalpy  change,  v^  (?  -  ?v)  ,  when  the  pressure  decreases  in 
the  exit  stream  to  the  saturation  value.  Altogether,  then, 
the  exit  velocity,  /  2  ih  +  g  az',  is 


■  {■  [■ 


C  (t  -  m  \  _  c  T  in  —  -1-  (?  -  o  )  /a 

PL  k  *s;  PL  s  T  ' rT  V;/3L 


*9<2I.-ZLh> 


']}l/2 


(A. 12) 


The  mass  flow  expression,  equation  (II. 7a)  follows  directly. 

The  density  of  the  exit  mixture,  for  an  isentropic  dis¬ 
charge,  would  be  similarly: 


^  -1 

{VLS  *  VpL  (tV,  -  VLS>  «  ^As}'1 


(A. 13) 


However,  previous  tests  (e.g.,  [25])  have  shown  that  the  evap¬ 

oration  from  the  liquid  occurs  mainly  at  the  surface  of  the 
stream  and  is  a  function  of  the  length-to-diameter  ratio  of 
the  orifice.  The  net  effect  is  that  the  exit  mixture  density 
has  an  average  value  that  is  larger  than  that  predicted  by 
equation  (A. 13) .  For  that  reason,  the  empirical  factor  3  has 
been  inserted  in  equation  (II. 7b).  Presumably,  a  factor  would 
also  be  needed  ..  the  velocity  relation,  equation  (A. 12),  but 
the  results  of  [25]  indicate  otherwise. 
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The  results  for  a  vapor  discharge  can  be  derived  simi¬ 
larly,  but  for  a  vapor  the  entropy  change  sv  -  sVs  is 


-  s 


Vs 


Z  R  in 


(A. 14) 


Equation  (II. 8a)  then  follows. 

A. 3  Evaporation  Rate  Correlation 

3y  observing  the  trends  exhibited  by  the  tests,  it  seemed 
that  the  larger  the  evaporation  rats  was,  the  larger  the  dis¬ 
crepancy  was  between  the  measured  tank  pressure  and  the  liquid 
saturation  pressure.  Consequently,  a  model  of  the  form: 


»W 


(we’sat  -  <we»test 


<we> 


n 

sat 


(A. 15) 


was  assumed.  (We) gat  could  be  calculated  for  any  test  from  the 
measured  liquid  temperatures  and  vapor  space  volumes;  (We) test 
could  also  be  computed  from  the  measured  vapor  space  pressures, 
temperatures,  and  volumes.  The  following  table  shows  some 
typical  numerical  results  (in  g/sec) : 


Test 

No. 

(We>sat 

(We} test 

AW 

e 

AW 

e 

'“e’sat 

AW  ' 
e 

1W<4at 

AW 

6 

2.43 

2.31 

0.12 

0.050 

0.020 

0.039 

20 

3.27 

3.17 

0.10 

0.031 

0.009 

0.023 

22 

0.079 

0.073 

0.0011 

0.014 

0.176 

0.026 

The  value  n  =  1.25  best  correlates  the  wide  range  of  evapora¬ 
tion  rates  shown  in  the  table.  To  make  equation  (A. 15)  dimen¬ 
sionally  correct,  the  right  hand  side  must  be  divided  by 
{ A_p T  /  RT'}  ^/'4  .  The  proportionality  constant  was  found  from 
linear  regression  to  be  1.9. 


A. 4  Derivation  of  Water  Ingestion  Model  Equations 

Derivation  of  equation  (II. -Qb)  .  When  the  relief  valve 
is  j armed  shut,  the  vapor  space  gas  is  compressed  by  the  in¬ 
gested  water  floating  on  the  cargo;  see  Figure  II. 2c.  Thus; 


(A. 16) 


where  y  =1  for  isothermal  compression,  and  Y  =  ratio  of  gas 
m  m 

specific  heats  for  adiabatic  compression. 

The  pressure  differential  at  the  puncture  is: 


IP 


»»  <z’- 


WT 


ZL) 


+  50L  (Zl  ~  3^) 


(A. 17) 


Since  it  is  required  that  d(AP)/dt  =  0,  equations  (A. 16)  and 
(A. 17)  can  be  combined  to  give: 


[  ?T  \  dZW 
\ZT  -  zwt;  dt 


0  (A. 13) 


But  dZm/dt  =  (Ww/ow  -  WLq/ol)/At  and  dZL/dt  = 

Equation  (11.20)  is  derived  by  substituting  these  relations 
into  equation  (A. 13)  and  solving  for  Ww> 

Derivation  of  equation  (11.23).  When  Z._m  <  Z„,  the 

‘  X  t  —4  • 

cargo  outflow  is  assumed  to  begin  when  the  water  level  in  the 
tank  is  slightly  greater  than  ;  in  particular,  it  begins 

at 


2., 


WT 


=  1.1  (Z. 


WT ,Eq . 


-  2  ) 
Lh 


+  Z 


Lh 


(A. 19) 


This  is  analogous  to  the  use  of  ?_.  -  ?_  _  and  Z_ .  -  Z.  _ 

XX  X  f  Zd ^  •  XjX  u  f  Lu 

in  equation  (11.17),  but  here  the  A?  at  the  puncture  is  just 


ZWT  "  ZLh  "  (ZWT,Eq.  *  ZLh)  0,1  (ZWT,Eq. 
(11.23)  ^follows  directly. 


ZLh> 


Equation 
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A. 5  Derivation  of  Choked  Flow  Equations  (11.31) 

As  shown  previously  in  equations  (A- 9)  and  (A. 10) ,  the 

thermodynamic  quality  of  the  discharge  mixture  for  an  isen- 

trooic  expansion  to  teiroerature  T  is: 

c 

*  *  CPLT=InF'/'c  <A'20> 

c 

a  T 

Making  the  substitution  T  =  T  (1  -  )  /  and  expanding  in  powers 

of  aT/T  gives: 

-  CPLiT  T1-  I  (t)  +  •'•]  (A'2l) 

Likewise 

“  ^Lc  =  CplaT  (A. 22) 

The  vapor  density  vVc  can  be  written  as: 


v 


Vc 


'V 


1  + 


AP 

P 


(A. 23) 


From  the  Chemical  Properties  File,  the  pressure  P  is  given  by 


^  £  =  A  —  — 

m  lo  A  (c  +  T") 


so 


AP 

P 


in  10 


B  AT 


(C  +  T*) 


There fore , 


the  product  *cvVc 


(A. 24) 


(A. 25) 


i 


can  be  expressed  as 


6 


X=VVC 


C  '  T 
^PL  1 


1  - 1  (¥)]  -'v  { 


AT 


1  +  TT 


3T  in  10 
(C  +  T*) 2 


i  ! 

J 


CPLiT  VV 


f ,  ^  IT 

3T  in  10 

3_ 

V  ~ 

_(C  +  T*) 2 

2 

(A. 26) 


Since  3xc  is  generally  snail,  the  tern  v^c  -  in  the  de¬ 

nominator  of  equation  (11.30)  is  approximated  as  v.  =  v^. 

Likewise,  X  *  X  ,  so  altogether  ecuation  (11.30)  can  be  writ- 
c 

ten  as : 


v 


L 


I  CPL 


M 


3CplV 


m  [> 


(A. 27) 


where  E  and  [D]  are  given  by  equations  (II. 31b)  and  (II. 31c). 
Equation  (A. 27)  is  quadratic  in  UT/T)  and  can  be  solved  to 
give  equation  (II. 31a). 


APPENDIX  B.  CHEMICAL  PROPERTIES  FILE 

Units  are  generally  in  the  cgs  system:  pressure  in  mul¬ 
tiples  of  105  dyne/cm2  =  102  N/m2  =  1  KPa ;  specific  heat  in 
cal/gr-°C;  latent  heat  in  cal/gr;  temperature  in  °C;  and  de.nsit 
in  gr/cm2.  The  temperature  range  of  validity  for  the  correla¬ 
tions  generally  correspond  to  vapor  pressures  from  about  0.1  to 
2  to  4  atmospheres.  Maximum  errors  of  the  correlations  with 
respect  to  the  cited  data  sources  are  also  given.  Correlations 
from  the  present  CPF  are  given  in  parentheses,  as  well  as  their 
maximum  error  over  the  same  temperature  range. 

In  most  of  the  data  sources,  English  units  are  used.  The 
new  correlations  were  initially  developed  therefore  in  English 
units  and  then  converted  to  cgs.  For  the  same  reason,  the 
English  units  correlations  given  in  the  CPF  were  reviewed  and 
converted  to  cgs. 

1.  ANHYDROUS  AMMONIA  [16,17,13],  molecular  weight  =  17.698 

Range  of  validity:  -40°C  to  38°C  (71.8  kPa  to  1461  kPa) 

97  6  98 

•  lQg10Pv  =  5  * 617 24  "  T  +  24  5.19  ?  nax*  error  3  -0.73% 

(=  7.1465  -  t^27  3~"2~  '  max'  error  =  4.05%) 

•  3X,  =  0.6371  -  1.413  X  10  2T  ;  max.  error  =  -0.29% 

(=  0.6533  -  1.299  X  10_2T  -  1.36  X  10-5T2  ;  max.  error  = 
2.43%) 

•  Z  =  0.932  -  1.647  X  10~3T  -  1.215  X  10_5T2  ;  max.  error  =  0 

(=  1.0  ;  max.  error  =  17.4%) 

•  C_  =  1.104  +  1.20  X  10  2T  ?  max.  error  =  0.73% 

t  la 

(1.100  +  1.000  X  10  2T  ;  max.  error  =  1.56%, 


3 


- 5  2 

.  c  =  0.257  -  4.02  X  10  T  -  1.504  X  10  T  ;  max .  error  = 

?v 

-2.61%) 

(Cpv  =  0.533  +  2.616  X  10~3T  +  3.756  X  10~7T2 
+  1.644  X  lO-13  ;  max.  error  =  449.9%) 

0.53741  P  (T  +  273.2) 

_  - - — - - =—  ,*  max.  error  =  -0.31% 

(T  +  24  5.2)2(0.6  371-1.413  X  lO^T) 

(=  327.0  ;  max.  error  =  13.82%) 

2.  BUTANE  [15,16,13,19,20] ,  molecular  weight  =  53.121 

Range  of  validity:  -55°C  to  43°C  (6.9  kPa  to  415  kPa) 

•  lo^10PV  =  5'9496  '  tI%|67'  ;  maX-  err°r  *  °'49% 

(=  5.9545  -  ;  max.  error  =  0.76% 

.  dt  =  0.5952  -  9.267  X  10~4T  ;  max.  error  =  0.14% 

L 

(=  0.6009  -  9.000  X  10_4T  ;  max.  error  =  1.37%) 

v  Z  -  0.95746  -  1.1405  X  10-3T  -  1.0404  X  lo'^T2  ;  max.  error 
*  -0.38% 

(=  1.0  ;  max.  error  =  12.55%) 

•  =  0.5540  +  2.7378  X  10-4T  +  1.5116  X  10_5T2  ;  max.  error 

•  PL 

=  1.64% 

(=  0.5800  +  1.300  X  10-3T  ;  max.  error  =  3.02%) 

•  =  0.3  35  ;  max.  error  =  3.7  2% 

(Cpv  =  0.3811  +  1.150  X  10“3T  -  6.384  X  10~7T2 
+  1.445  X  10-1°T3  ;  max.  error  =  23.38%) 


X  =  [7 


1.069  -  8.4656  X  10"2T  -  7.7226  X  10~4T2  j  3-y ) 


0.51897  Pv  (T  +  273.2) 


(T + 239.7) 2(0. 5952  -  9.267  X  10~4T 
(=  92.0  ;  max.  error  =  13.58%) 


;  max.  error  =  1.57% 


3.  CHLORINE  [21],  molecular  weight  =  70.91 


Range  of  validity:  -51°  to  38°  (45  kPa  to  1067  kPa) 

.  log1QPv  =  6.054  0  -  *  max.  error  =  0.16% 

(=  6.5425  -  f j'1 1 >  max*  error  =  5.02%) 

«  =  1.4676  -  2.8738  X  10  3T  ;  max.  error  =  0.49% 

(=  1.4599  -  2.599  X  10~3T  ;  max.  error  =  0.93%) 

•  Z  =  0.948  -  1.3064  X  10-3T  -  9.383  X  10_6T2;max.  error 
=  0.10% 

(  =  1.0  ;  max.  error  =  11.5%) 

•  C?L  =  0.2334  +  1.35  X  10  4T  ;  max.  error  =  -2.05% 

(=  0.1200  +  5.00  X  1Q_4T  ;  max.  error  =  -58.19%) 

•  C77I  =  0.0332  -  5.6571  X  10_4T  -  4.3393  X  10_6T2  ;  max.  error 


PV 
=  0.61% 

(Cpv  =  0.1130  +  5.8933  X  10_5T  -  9.166  X  10~8T2  ; 
max.  error  =  89.73%) 


=  [52.585  -  7.2465  X  10_2T  -  5.2047  X  10~4T2jf 

J\ T  +  24  6.02/ 


0.47317  Pv  (T +  273.2) 


(T  +  2  46. 52)  2  (1.4  67  6  -  2.87  38  X  10-2T) 
(=  63.72  ;  max.  error  =  19.28%) 


;  max.  error 


-0.55% 


4.  LNG  ii5,13,19,22(23],  molecular  weicht  =  16. Cl 


Range  of  validity:  -177°  to  -153°  (22.3  kPa  to  202.5 


.  „  ,  389.39 

1  g10‘ V  5.736a  -  T  +  266 . 0 

( =  same ) 


;  max.  error  =  -1.14% 


=  0.1333  -  1.665  X  10  T  ;  max.  error  =  0% 

(=  0.1880  -  1.399  X  10_3T  ;  max.  error  =  1.98%) 

Z  =  -2.5724  -  3.97015  X  10~2T  -  1.10412  X  10~4T2  ;  max.  error 
=  0.45%) 

(=  1.0  ;  max.  error  =  4.21%) 


=  1.1995 

4* 

2.2639 

X 

10~3T  ; 

max.  error  = 

-0.57% 

(=  1.1961 

+ 

2.1996 

X 

10“3T  ; 

max.  error  = 

1.33%) 

=  2.2639 

+ 

1.0904 

X 

10“2T  ; 

max.  error  = 

-4.48% 

(CpV  =  0.5166  +  6.6838  X  10-4T  +  5.311  X  10~7T2;nax. 
error  =  12.1%) 


=  -  [236.63  +  4.4238  T  +  1.23023  X  10_2T2]  ( 2  ^ 


0.21447  P  (T  +  273 . 2) 


(T  +  266.0)  (0.1833  -  1.665  X10  JT) 

{=  121.8  ;  max.  error  =  -4.65%) 


;  max.  error  =  1.51% 


Assumed  to  be  methane. 
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METHYL  ALCOHOL  [2,16,13,19,24] 


molecular  weight 


32.043 


Range  of  validity:  -10°  to  80°  (2  k?a  to  179  kPa) 


log 


=  7.20535  - 


1582.27 


10  V 
(=  7.9265  - 


T  +  23  9.73  ; 

2002.22 


T  +  273.17 


max . 


max.  error  =  -5.45%  at  -10° 
max.  error  =  0.32%,  T  >  15° 

error  =  2.86%) 


=  0.8099  -  9.1979  X  10  4T  ;  max.  error  =  0.10% 
(same) 


Z  =  0.959  +  5.000  X  10~4T  -  1.0833  X  10“5T2  ;  max.  error  = 
=  0.13%) 


(= 

1.0  ;  i 

max.  err 

or 

=  7.53%) 

II 

■4 

Cl. 

u 

0.5660 

+  1.622 

X 

10“3T  ; 

max . 

error  =  -0.43% 

(= 

0.5260 

+  3.200 

X 

10_3T  ; 

max . 

error  =  9.26%) 

0.3018 

-  1.601 

X 

1 

m 

l 

o 

rH 

1.905 

X  10  V  ;  max.  error 

=  3.00% 

(Cpv  =  0.2922  +  6.886  X  10"4T  +  2.497  X  10_7T2  ;  tax. 
error  =  100.72%) 


X  =  [216.56  +  1.129  X  10_1T  -  2.446  X  10  3T2 


T  +  273 . 2 y 


0.87037  Pv  (T  +  273.2) 

(T +  239.7) 2  (0. 8099  -  9.1979  X  10~4T) 


T  +  2  3  9  .  /  / 

;  max.  error  =-1.96% 


(=  262.8  ;  max.  error  =  -6.58%) 


PROPANE  [15 


,19,20,25  ],  moleculai 


»e:cr.t  =  3  3 


Range  of  validity:  -34°  to 


lo<?10?v  =  5-9633  -  tT 


3  IS . 56 


243.63  ; 


(=  5.9545  - 


313.33 
T  +  2  4  3  .  0  ; 


max . 


44°  (1.5  kPa  to  1520  k? 

max.  error  =  -1.67%  at 
max.  error  =  -0.73?,  T 


error  =  -1.20%,  T  <  27° 


cL  =  0.5301  -  1.425  X  10_3T  -  2.738  X  10~6T2  ;  max.  err 
=  0.46% 

(=  0.5350  -  1.100  X  10  3T  ;  max.  error  =  5.65%) 


2  =  0.8395  -  2.499  X  10-3T  -  1.436  X  1Q~5T2  ;  max.  erro 
=  0.85% 

(=  1.0  ;  max.  error  =  23.3%) 


C?L  =  0.5731 
=  -6.46% 

(=  0.6599 


1.6913  X  10  ■'T  +  1.7  05  X  10  3T2;  max.  err 


+  2.50  X  10  3T  ;  max.  error  =  -15.41%) 


-  1.6045  X  10~3T  -  1.4134  X  10_5T2; 

-21.23%  at  -70° 

2.52%,  T  >  -70° 

(C„7  =  0.3700  1.153  X  10"3T  -  4.083  X  10~7T2;  max 

r  V 

error  =  154.32%) 


-PV  -  °-  2373 
max.  error  = 
max.  error  = 


X  =  [75.512  -  2.121  X  10_1T  -  1.2615  X  10_3T2] 
0.45026  ?v  (T  4-  273.2) 

-  - m - —  ;  rax .  e: 

(T  +  248.7)  (0.5301  -  1.10  X  10_JT) 


/  T  +  27  ; 
VT  ”  24: 


:or  =  ^ . 


(=  101.78  ;  max.  error 


46.72%) 
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APPENDIX  C.  COMPUTER  PROGRAM 


The  analytical  model  presented  in  Section  II  has  been 
programmed  for  numerical  solution  in  the  computer  code  TVENT, 
which  is  outlined  in  this  appendix.  The  basic  technique  used 
in  the  code  is  to  assume  that  a  small  percentage  of  the  cargo 
mass  is  discharged  and  then  compute  the  time  increment  for  the 
discharge  to  occur,  using  the  various  energy  and  mass  balances 
and  flow  relations.  After  the  model  has  been  "solved"  to  a 
prescribed  degree  of  accuracy  for  this  incremental  discharge, 
am  additional  discharge  is  assumed,  and  the  solution  process 
repeated.  This  continues  until  tank  conditions  are  computed 
that  make  further  discharges  impossible. 

TVENT  is  arranged  into  a  "MAIN"  code  and  a  number  of  sub¬ 
routines  that  generally  correspond  to  the  various  submodels 
comprising  the  entire  venting  rate  model.  Flow  charts  giving 
the  logic  for  the  MAIN  program  and  the  subroutines  are  given 
in  Figures  C.l  through  C.12  and  Tables  C.2,  C.3,  and  C.4. 

Table  C.l  gives  the  required  inputs  for  the  program  in  its 
present  form. 

To  describe  the  code  briefly,  the  MAIN  program:  (1)  or¬ 
ganizes  the  input  data;  (2)  determines  whether  the  cargo  is 
volatile;  (3)  determines  whether  liquid  or  vapor  or  both  are 
discharged;  (4)  selects  the  increment  of  mass  to  be  discharged; 
(5)  determines  whether  the  calculations  for  any  one  iteration 
have  converged;  (6)  determines  when  the  discharge  is  completed; 
and  (7)  calls  the  various  subroutines  that  perform  the  detailed 
computations.  Before  each  calculation  cycle,  MAIN  also  calls 
the  subroutine  INGEST,  which  determines  whether  the  tank  con¬ 
ditions  are  appropriate  for  air  or  water  ingestion.  A  test 
for  choked  flow  conditions  (CHOKTEST)  is  called  before  the 
first  cycle  and  each  cycle  thereafter  so  long  as  choked  out¬ 
flow  can  occur. 
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Subroutines  are  used  to  make  most  of  the  calculations. 
ENERGY  solves  the  energy  and  mass  balances  for  a  volatile 
cargo  and  calls  FL0W1  to  compute  discharge  flow  rates.  NOMENER 
and  FL0W2  are  used  similarly  for  nonvolatile  cargos.  INGEST 
performs  all  the  logic  operations  needed  to  determine  whether 
air  or  water  ingestion  can  occur  at  any  instant.  AIRIN  com¬ 
putes  the  change  in  tank  conditions  when  air  bubbles  are  in¬ 
gested  by  calling  ENERGY2  (for  volatile  cargos)  or  N0NENER2 
(for  nonvolatile  cargos)  to  determine  the  vapor  space  pressure 
rise;  a  return  to  ENERGY  or  NONENER  is  then  made  to  compute 
the  subsequent  liquid  discharge.  WATIN  computes  the  cargo 
discharge  and  water  inflow  when  water  is  ingested,  and  calls 
FL0W2  to  compute  the  outflow.  NOOCIT  is  used  to  compute  the 
discharge  when  water  can  be  ingested  at  the  instant  the  punc¬ 
ture  is  opened,  as  well  as  for  certain  other  water  ingestion 
cases  similar  to  this.  Other  self-explanatory  subroutines 
compute  cargo  properties  and  miscellaneous  data. 

A  computer  listing  of  each  routine  is  also  given  in  this 
appendix.  Table  C.5  explains  the  correspondence  between  the 
symbols  used  in  the  analysis  and  those  used  in  the  computer 
code.  Table  C.6  gives  the  formats  to  be  used  to  enter  the  in¬ 
put  data  in  TVENT. 

Other  notes  on  the  use  of  the  computer  code  include: 

1.  The  normal  outputs  are  the  total  mass  of  cargo 
liquid,  cargo  vapor,  and  air  released,  and  the 
total  discharge  time. 

2.  The  code  assumes  that  only  liquid  or  only  vapor 
can  be  discharged  through  circular  or  horizontal 
slot-like  punctures;  when  the  liquid  surface  inter¬ 
sects  this  type  of  puncture,  only  liquid  is  dis¬ 
charged  until  the  level  drops  y  -ow  the  puncture. 
Simultaneous  liquid  and  vapor  discharges  are  allowed 
for  vertical  slot-like  punctu  es  when  the  liquid 
level  intersects  the  puncture 
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3.  Similarly,  circular  or  horizontal  slot-like  punc¬ 
tures  must  be  either  completely  open  to  the  air  or 
completely  submerged.  Only  for  a  vertical  slot¬ 
like  puncture  is  the  water  level  allowed  to  inter¬ 
sect  the  puncture. 

4.  The  code  assumes  the  cross-sectional  area  of  the 
tank  is  constant;  for  spherical  or  other  tank 
shapes,  an  effective  tank  height  should  be  input 
such  that  Vip/ZT  gives  a  representative  or  average 
area. 

5.  For  nonvolatile  substances , only  the  liquid  density 
correlation  must  be  input  accurately.  The  other 
correlations,  if  they  are  not  known  for  that  cargo, 
can  be  input  as  any  convenient  values  (say,  one) . 
Care  should  be  taken,  however,  to  make  sure  that 
the  fictitious  vapor  pressure  correlation  predicts 
positive  pressures  that  are  significantly  less  than 
atmospheric. 
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COMMON  au's 

SYMBOL  DEF'S 
AND 

CONSTANT 

TYPE  DBF'S 

OPTION 

DESCRIPTIONS 

DEF'S 

I  CARGO,  PUNC¬ 
TURE,  SAFETY 
VALVE  INPUT 
AND  AIR 

TEMP,  DEPIN1TI0N 


CHEMICAL 

PROPERTIES 

INPUT 


WATER  HEIGHT 
AND 

SPECIFIC 

VOLUME 


SEE 
FLOW 
CHART 
TVENT,  A 
FOR  ALGOR 
DETAILS 


♦ 

230 

INITIAL  CONDI¬ 
TIONS,  INITIAL 
PRESSURES 
DETERMINATION 

CONTINUE 

INITIAL 
PARAMETERS 
DETERMINATION 
AND  SPECIFICA¬ 
TION 

INITIATION 

AND 

CALL  INCEST 
DHL  AND  DMV 
SPECIFIED 

DHL  AND 
DMV  ARE 
LIQUID 


INCREMENTS , 
RESPECTIVELY 


INGESTION 


(SEE  TVENT.  B 
FLOW  CHART) 


FINAL  OUTPUT 
AND  FORMATS 
FOR  TVENT 


CARGO 

RELEASE 

LOOP 

Ui&RliI 

AND  FLOW 

EQUATIONS 

I.OOP 

DETAILED  FLOW  CHART  DIAGRAMS 
FOR  THE  CARGO  RELEASE  LOOP 
AND  FOR  THE  ENERCY  AND 
FLOW  EQUATIONS  LOOP 
ARE  GIVEN  ON  FLOW  CHARTS 
TVENT,  B  AND  TVENT,  C, 
RESPECTIVELY. 


PARAMETER  SUBSCRIPTS 


6i  +  4t  or  Ei+1 


FIGURE  C-l.  PROGRAM  TVENT  (MAIM  PROGRAM)  FLOW  CHART 


(Ill, 


600 


FIGURE  C-I.  PROGRAM  TVENT  (MAIN  PROGRAM)  FLOW  CHART  (Cont’d) 
(TV ENT,  3  -  Cargo  Release  Loop  Flow  Chart) 


FIGURE  C-2.  ENERGY  SUBROUTINE  FLOW  CHART 


W  -  CpL(T  -  Tj.)  -  TS.(S  -  SL.)  +  104(Pt  -  Pv)VI>/J  +  -  ZUi)/J 

**DUM1  -  2(104(Pt-PJVl  +  4(2^  -  Z^)) 


FIGURE  C-4.  SUBROUTINE  FLOW1  FLOW  CHART 


CALCULATE 
KLCi),  ZU4) 
ZLB 


107 


600 


FIGURE  C -8.  SUBROUTINE  CHOKTST  (KTYP)  FLOW  CHART 
(Concl 'd) 


970 


FIGURE  C-9.  SUBROUTINE  INGEST  FLOW  CHART  (Cont'd) 


22 


TEST  -  PATH  ♦ 
10-4C(ZW-ZL(4))/VW/2 


ICT  -  ICT  -  1 
SET  TIME  STEP,«t 


2213 


CALCULATE 

WLO 


ML (4)  -  ML (2)  - 
WLO*DT 

SCROLL  PARAMETERS 
CALCULATE  TOTAL 
MASS  LOST 


FIGURE  C-9.  SUBROUTINE  INGEST  FLOW  CHART 


TEST  -  PATH  ♦ 

lO^CCZV-ZP+DP)/ 

'Till 


ICT  -  ICT  -  1 
CALCULATE  CAS 
PARAMETERS 


~T~ 

CALCULATE 
TIME  STEP 


2311  f  ' 

CALCULATE  CAS 

THERMOPHTSICAL 

PARAMETERS 

V 

CALCULATE 
GAS  PLOW 
RAXES 

■  j 

CALCULATE 
TOTAL  MASS 
LOST 


SCROLL  PARAMETERS 


FIGURE  C-IO.  SUBROUTINE  AIRIN  (ITIME)  FLOW  CHAR 
(Concl'd) 


SVL(4)  >  VW 


CALCULATE  INITIAL 


VALUES  OF  PTINT 


ZLINT,  FT,  T,  ZVT 


CALCULATE  FLOW 


AREAS  AND  FLOW 


SATE  (VW-WLO) 


CALCULATE  DT 


ML(4),  ZL(4) 


ZLH,  ZVT(4) 


|WW-OLDWW| 


>  c • OLDWW 


AND  ACCUMULATE 


ZL(4)  >  ZP-DP 


121 


500 


FIGURE  C-ll.  SUBROUTINE  WAT IN  FLOW  CHART 
(Cone 'd) 


i 


300 


CALCULATE 

INITIAL  VALUES 

OF  PTINT,  ZLINT, 

PT,  ZVT,  ZL 

L  _ 

r 

CALCULATE  PT(3)  •  ?T(4), 
WLO ,  WW,  DT,  ML(4), 
ZL(4),  ZLa,  ZWT(4), 
PT(4),  T(4) 

KM 

Mi 

SCROLL  PARAMETERS 

AND  ACCUMULATE 

LOST  MASS 

FIGURE  C-ll.  SUBROUTINE  WAT IN  FLOW  CHAR' 
(Concl’d) 


L23 


FIGURE  C-12.  SUBROUTINE  NOOUT  FLOW  CHART 


FIGURE  C-12 .  SUBROUTINE  NOOUT  FLOW  CHART 
(Cont 'd) 


FIGURE  C -12 .  SUBROUTINE  NOOUT  FLOW  CHART 
(Cont  Td) 
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TABLE  C.l.  ANTING  RATE  MODEL  INPUT 


1. 


3. 


5. 


6. 


7. 


8. 


9. 


10. 


11. 


12. 


13. 


14. 


15. 


16. 


17. 


18. 


Cargo  Identlf ication  code  (to  be  used  eventually  with  the  HACS 
Chemical  Properties  File) 

Cargo  molecular  weight 

Total  initial  mass  of  cargo,  grams 

Cargo  initial  temperature,  °C 

Volatility  constant;  used  in  determining  whether  the  energy  balance 
for  volatile  or  nonvolatile  cargos  will  be  used,  by  comparing  satura¬ 
tion  pressure  at  the  cargo  initial  temperature  to  (volatility  constant) 
times  (atmospheric  pressure);  a  value  of  0.8  has  been  found  acceptable. 

Air  temperature,  aC 

Air  pressure,  KPa  (Standard  Atmosphere  »  101.4  KPa) 

Tank  volume,  cm^ 

Tank  height,  cm 

Estimate  of  tank  pressure,  KPa;  is  compared  to  atmospheric  pressure 
and  initial  estimate  of  tank  pressure  is  set  to  larger  of  the  two. 

Puncture  description:  round,  horizontal,  or  vertical  slot 

Height  of  puncture  above  tank  bottom,  cm 
2 

Puncture  area,  cm 

Puncture  width,  if  horizontal  or  vertical  slot,  cm 

Puncture  discharge  coefficient 

Safety  valve  condition:  stuck  or  unstuck 

Safety  valve  vacuum  differential  setting,  if  unstuck,  KPa 

Height  of  outside  water  above  tank  bottom,  cm 


19-24.  Constants  in  thermodynamic  property  correlations  for  vapor  pressure 
(KPa),  liquid  density  (g/cm^),  vapor  compressibility  (-) ,  liquid 
specific  heat  (cal/g-’C)  ,  vapor  specific  heat  (cal/g-*C) ,  and 
latent  heat  of  evaporation  (cal/g). 

25.  Descriptor  for  adiabatic  or  isothermal  conditions. 
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TABLE  C-2 .  MISCELLANEOUS  SUBPROGRAMS 


SUBROUTINE  AISMASS : 

Uses  ideal  gas  equation  of  state  to  determine  mass  of  air  in  the 
cargo  tank  vapor  space. 


SUBROUTINE  CORRECT: 

Updates  latent  heat,  specific  volume  of  vapor,  specific  volume  of 
liquid,  liquid  specific  heat,  and  vapor  specific  heat  with  the  most 
recent  value  of  temperature  at  time  t^+^. 


FUNCTION  CPLHF  (T) : 

Computes  the  latent  heat  of  vaporization  of  the  cargo  liquid  with 
temperature,  T. 


FUNCTION  GAMH  (RA) : 

Computes  the  average  value  of  the  ratio  of  specific  heats  of  the 
vapor  space  gases  (cargo  vapor  and,  when  present,  air).  The  average 
is  taken  over  times  ^  and  t^. 


SUBROUTINE  TERMIN  (END  CODE)  : 

A  terminal  printout  of  time,  mass,  pressure,  height,  and  temperature. 
TERMIN  is  entered  from  another  subprogram  that  is  identified  by  the 
name  given  by  the  Hollerith  Field  designated  by  the  parameter  end 
code. 


SUBROUTINE  ZLHCALC : 

Computes  the  flow  center  height,  ZLH. 
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TABLE  C-3  .  CPF  ( IFLAG, T)  FUNCTION  SUBPROGRAM 


IFLAG  -  Inceger  co  describe  which  thermophysical  property  is  to  be 
determined 

T  -  Independent  parameter 


IFLAG 

T 

CPF 

1 

Cargo  Temperature 

Cargo  Material 

Vapor  Pressure,  P  ■  f(T) 

2 

Cargo  Temperature 

Cargo  Material 

Liquid  Density 

3 

Cargo  Temperature 

Cargo  Vapor  Compressibility 
Factor,  Z 

4 

Cargo  Temperature 

Cargo  Material 

Liquid  Specific  Heat 

5 

Cargo  Temperature 

Cargo  Material  Saturated 
Vapor  Specific  Heat  at 
Constant  Pressure 

6 

Cargo  Material 

Vapor  Pressure 

Cargo  Temperature, 

T  -  f_1(P) 
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TABLE  C-4.  SUBROUTINE  INGEST  FLAG  PARAMETERS 


One  key  parameter  that  determines  the  computational  process  is  called 
INGFLG.  Values  that  INGFLG  can  have  and  the  effect  on  computation  pro¬ 
gression  are  given  below: 


INGFLG  RESULT 


0  Uninterrupted  venting  -  normal  mode.  Adequate  pressure 

within  the  tank  for  liquid  cargo  and/or  gas  (cargo  vapor 
and  possibly  air)  to  vent  through  the  tank  puncture. 

2  Insufficient  pressure  within  the  tank  to  vent  to  the  out¬ 
side  media  (air  and/or  water) .  Initiated  in  subroutine 
FL0W1  or  FL0W2.  Used  as  a  transfer  flag  to  route  computa¬ 
tion  process  to  subroutine  INGEST.  INGFLG  will  be  changed 
to  another  appropriate  value  in  INGEST. 

1  Initiated  in  subroutine  INGEST  to  route  computation  process 

to  subroutine  AIRIN.  If  INGFLG  is  set  to  1  in  INGEST,  return 
from  INGEST  to  the  main  program,  TVENT,  is  followed  by  a  call 
to  subroutine  AIRIN.  Return  to  TVENT  from  AIRIN  causes 
INGFLG  to  be  reset  to  0. 

3  Initiated  in  subroutine  INGEST.  Causes  return  to  the  main 
program,  TVENT,  and  an  immediate  call  to  subroutine  WATIN. 

4  Initiated  in  subroutine  INGEST.  Causes  return  to  the  main 
program,  TVENT,  and  an  immediate  call  to  subroutine  NOOUT. 

Initiated  in  subroutine  INGEST.  Causes  air  to  be  admitted 
to  the  vapor  space  in  the  tank  instantaneously.  Return  to 
the  main  program,  TVENT,  follows,  whereupon  in  TVENT,  INGFLG 
is  reset  to  0  and  the  cargo  release  loop  is  reentered. 


5 
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TABLE  C.5.  LIST  OF  VARIABLES  IS  COMPUTER  PROCRAM 


7ariabla 

AISI 

Cal-» 

I*b«l 

INC 

Ocscrlsclan 

Saucia«a  Vh*r«  3«f«r«nc*<i 

Air  bubble  volua* 

«•*£  ict*nc 

were,  add 

AO C 

•* 

ca* 

G2* 

Araa  of  jaa  flew 

west,  new,  now, 

HOIST 

AOL 

2 

cm 

GO 

Araa  of  liquid  flow 

west,  now,  now 

A? 

INC 

Araa  of  puneeura 

west,  add,  soout. 

•watd 

at 

2 

ca 

CO 

Tank  croaa-aaceional 

araa 

WEST,  ADMASS ,  ESEKCT. 

herot:.  hoist,  scseser, 

S0KESC2,  SOOUT,  WATTS 

Cl.  Cl 

— 

C? 

Conacaaca  uaad  in 
calculaciag  apaclfle 
haac  of  cargo 

WDT.  C?T 

C301P  (C3) 

- 

CO 

Siacharga  eoaffi- 
ciaac 

were,  now.  now, 

HCEST ,  SOOUT,  ’WAITS 

Oil 

S 

GO 

laicial  cargo  saaa 

WEST 

OWT  (MOLE) 

3 

g  -  aola 

GO 

Cargo  aolacular 
waighe 

WEST,  ADMASS .  CEOKTST , 
CORRECT,  ESERGT,  ESERCT2, 

now,  now.  gams,  dcest, 

SOSESER,  S0SZXC2 

CPA 

cal 

S-’C 

GO 

S  pacific  haac  of  air 

were,  BWERGT.  ESERCT2, 

now,  now,  cams,  hoist 

CPL 

cal 

T=*3 

GO 

Specific  hear  of 
cargo  liquid 

were,  ADIS,  CORRECT, 
OIERGT,  ESERGT2,  CAMS, 
HOIST,  WATD 

CTT&AR 

cal 

GO 

Specific  haac  of 
cargo  vapor 

WEST,  AIRIS,  CORRECT, 
HERGT,  ESERCT2,  CASS!, 
HCEST,  'WAITS 

91,  02 

— 

C? 

Cotucanc*  uaad  is 
calculaciag  liquid 
dona ley 

WEST,  CPF 

del  time 

i 

INC 

Svaporacad  aasa  par 
eiaa  aca? 

WEST,  ESERCT,  ESERGT2, 
SOOUT 

CMA 

t 

OG 

Maaa  of  air  ralaaaad 
par  eiaa  acap 

WOT,  ADIS ,  ESERGT, 
SOSESER 

30, 

I 

OG 

Maaa  of  liquid  ra¬ 
laaaad  par  eiaa 
•cap 

were,  ADD,  ESERGT, 
HCEST,  SOSESER,  SOOUT, 
WATD 

TABLE  C.5. 


1ST  OF  VARIABLES  IN  COMPUTER  PROGRAM  (Conc'd) 


Inles 

Conaon 

Label 

Inscription 

Routines  'Where  Referenced 

Variabl# 

3JW 

t 

oc 

Mas*  of  vapor  re- 
leaaed  par  tlaa  seep 

TVOT,  AIRIS,  EtERGT. 
ISGEST,  SO  MOO 

DP 

cm 

rsc 

Vertical  extant  si 
vertical  aloes 

TV  ENT,  ENERCY.  ISCEST. 
SONESES,  NOOCT .  WATTS, 
ZLHCA L£ 

DPV  OELPVAL) 

CP* 

GO 

Safari  valve  relief 
vacuua  setting 

tvot.  oerct.  incest, 
SONENER 

orsub 

— 

INC 

Solaces  energy  auo- 
roue la* 

TVOT,  AIRIS 

rtosra 

— 

INC 

Solaces  flow  aub- 
roueln* 

tvot ,  airis,  grcetst 

a 

ca 

sac 

GO 

Gravitational 
accalarae ion 

rrart,  ma,  fxovx,  now? 
ISGEST,  SOOCT ,  WATIN 

31*97 

— 

C? 

Cooacancj  for  lacane 
haac  of  cargo 

rrot,  cplhp 

BID  (IHEATTR) 

— 

GO 

Haac  cransfar  da- 
scripeor 

rrot,  energy,  energy?, 
digest.  nosesex,  nonesc:, 

WATTS 

:ei 

_ 

rsc 

I la*  atop  counter 

TVEMT,  INCEST 

rscru; 

— 

OG 

Subroutine  flag  need 
in  IHCtST 

tvot,  axro.  oergy, 
FL0W1,  now?,  INCEST , 
SONESER 

s 

1 

■  «-<=»■  ■ 

GO 

Mechanical  equlva- 
lane  of  haac 

T7ENT,  AI3MASS,  FL0W1 

eal-aac* 

LAMBDA 

cal 

i 

CO 

Latent  haac  of  cargo 

TVOT.  CORRECT ,  ENE3CY, 
ESERZ 

.HA 

t 

GO 

Air  name  la  rank. 

TVOT.  AIRIS ,  A1  AMASS, 
ESERCY.  ENE3GYI,  FLOtfl, 
PLCW2.  GAJCt,  INGEST , 
NONINEK,  NON ENC2,  SCO LT. 

termo,  watts 

ML 

i 

GO 

Liauid  naee  In  rank 

TVOT,  AIRIS.  A  IRMAS  S  , 
CORRECT,  ESERCY ,  ESERCY?. 
ISGEST.  SONENER.  NONESC 2, 
sooerr,  tesmo,  «ato 

MV 

( 

GO 

7apor  aaaa  la  lank 

TVOT,  AIRIS,  AIRMASS, 
CORRECT,  ENERGY,  ESERCT?, 

CORRECT,  CTERGT,  OERCT?, 
flcwi,  n.ow2,  camk, 
ntcESx,  nonoer,  nonoci, 
xooot,  rsuns,  watts 
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TABLE  c. 

. G .  LIST 

of  var: 

Cowon 

A3LZS  IN  COMPUTER 

7*ri*bi* 

lairs 

Laos! 

Gsscripcion 

.ms? 

- 

INC 

Loop  couotsr  through 
aaia  progrra 

?i  -  ?3 

— 

C? 

Constants  i or  ralcu- 
Lacing  /spar  prsssurs 

?A 

CPs 

CDJ 

Air  prsssars  in  chs 
cans 

?AT2T 

CPs 

SEN 

Aoosphtrlc  prsssurs 

fZSST 

CPs 

SEN 

Aaoltcc  prsssurs  «x- 
csrnal  to  puncturs 

PINTS 

K?s 

GEN 

Aaolsnc  prsssurs  ix- 
csraai  to  puncturs 
for  gas  flow 

CPs 

GDI 

Tank  prsssurs 

pnqi 

K?s 

’•iA? 

Tank  prsssurs  sc 
squUlbrl-ja 

?v 

Os 

ca 

Vapor  prsssurs  la 
ebs  rank 

*c 

Os-ca3 

g-aols-'* 

ca 

Gnivsrsal  gas  con- 
scant  (  aschanlcsi 
units 

*c: 

Cal 

g-*ols-’X 

Ga 

Taivsrsai  gas  con¬ 
stant,  thsraal  units 

S7A 

3 

ca 

1 

ca 

Spsclfic  vcluas  of 
sir 

STD  (IVAL7I) 

— 

ca 

Safaey  t*J.vs  ds- 
scrlptor 

STL  (TL) 

si 

i 

ca 

Spsclfic  vol UBS  of 

Liquid 

PRXRA.M  (Cone ’a) 


RcutlstA  where  Referenced 
T.-ENT,  INCEST 


men,  choktst,  c?r 


r/orr,  airxn,  enesct, 

WAT  a 

TVEIT,  ENERCT,  INCEST, 
NCMENER 

ttsnt,  choktst,  r.owi, 
r.yn ,  nicest.  mocot 

TVENT,  CHOKTST,  51CV1 , 
TIC  K2 


merr,  airxn ,  aco-ass, 

CHOKTST,  CORRECT.  ENERCT, 
ENERGT2,  ELCUl,  FLOV2, 
ONCOST,  NCNTNER,  MCSENC2 . 
NOOCT,  TE3MIN,  WAT  IN 

ttent,  incest,  waits 


tvett,  a  ms,  axe-ass. 
CHOKTST,  CORRECT.  ENERCT, 
INERCT2,  INCEST ,  SCNENER. 
SONINC2,  NOOCT,  TERMEI, 
WATTS 

EVENT,  A  TEAS  S,  CHOKTST, 
CORRECT,  ENERCT,  -SE5CT:, 
nova,  flows,  CAXM.  INCEST, 
SOSESEH.  N0SENC2,  NOOCT, 
WATTS 

mm.  nowi,  now2,  cans-, 

INCEST 

ttest,  enerctt,  scsenc: 


T7ENT,  EMERGE,  NICEST, 
SONENES,  MOOOT ,  WATTS 

mm,  AXRIN ,  CORRECT, 
ENERCT,  ESERCT2,  INCEST, 
MOSESES,  MCNENC2,  NCCUT, 
WAITS 
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TABLE  C.5.  LIST  OF  VARIABLES  IN  COMPUTER  PROGRAM  (Cont'd) 


Variable 

Units 

Common 

Ub«l 

Description 

leucines  Where  Referenced 

SW  (VV) 

3 

ea 

i 

GE1 

Specific  volume  of 
vapor 

TVEMT,  AKIN,  CEOKTST, 
CORRECT,  EMERGE,  ENERGE2, 
INGEST,  WAIIM 

•c 

GIN 

Cargo  teoperacure 

TV  ENT,  AHUM,  AI3MASS, 
CHOKTST,  CORRECT,  CPF, 
CFLHF,  EMERGE,  EMERCT2 , 
FL0W1,  FL0W2,  GAMM,  INGEST 
NONENER,  NONENG2,  MOCCT, 
TESMTM,  WATIN 

TC 

•c 

GEM 

Initial  cargo  tem¬ 
perature 

T7ENT,  NOMEMSR,  N0NENG2 

TIME  (TIM) 

sec 

GEM 

Tlae 

TVEMT,  AIR  IN,  EMERGE , 
INGEST,  MOOUT,  TSSMIM, 
MATIN 

TIN?  (TATM) 

•c 

GEM 

Air  temperature 

TVEMT,  EMERGE,  EMERGT2 

IMA 

g 

IMG 

Total  aass  of  air 
released  from  tank 

TVEMT,  AKIN,  INGEST, 
MOOUT,  TERMIN,  WAT  IN 

TML 

t 

IMG 

Total  mass  of  liquid 
released  from  tank 

TVEMT,  AIRIN,  INCEST, 
NOOUT,  TERM  IN,  VATIN 

TMV 

t 

IMG 

Total  aass  of  vapor 
released  from  tank 

TVEMT,  AKIM,  INGEST, 
MOOUT,  TERMIN.  WAT  IN 

V1-V3 

~~ 

C? 

Constants  used  in 
calculating  specific 
beat  of  vapor 

TVEMT,  CPF 

7QU 

3 

cm 

IMG 

Volume  of  air  bubble 

TVEMT,  AIRIN,  ENERGT2, 
SOMENG 

VT  (VOLT) 

3 

ea 

GEM 

Tank  volume 

TVEMT ,  AKMASS ,  ENERCE, 
EHERCT2,  INCEST,  NONENER, 
MONENC2,  MCOCT,  VATIN 

TO 

3 

ea 

i 

GZM 

Specific  volume  of 
water 

TVEMT,  INGEST,  NOOUT, 

WAITS 

«AO 

X 

s*c 

GEM 

Discharge  race  of  air 

released 

TVEMT,  AKIM,  EMERGE, 
•MERCT2,  FL0W1,  FL0W2, 
INGEST,  NONENER,  N0NENG2 

WLO 

X 

*«c 

GSM 

Discharge  rata  of 
liquid  released 

TVEMT,  AIRK,  EMERGT, 
•ME3GT2,  FL0W1,  FLCW2, 
INGEST,  NONENER,  S0SENG2, 
NOOUT,  VATIN 

VVO 

1 

amc 

GEM 

Discharge  rate  of 
vapor  released 

TVEMT,  AKIM,  EMERGE, 
EMERGE!,  FLOW1,  FL0W2, 
INGEST,  NONESER,  NONEMG2 

i 

i 


TABLE  C .5 


LIST  OF  VARIABLES  IN  COMPUTER  PROGRAM  (Concl'd) 


GooBOn 


VarlabU 

Uni  :  s 

Label 

Description 

Routines  Where  Referenced 

XL? 

ca 

ING 

Sloe  width 

tvest,  digest,  soour 

21 -S3 

C? 

Constants  used  for 
calculating  vapor 
compressibility 

77DIT,  C7F 

a 

ca 

GSS 

Haight  o f  liquid  in- 
slda  tank 

T7E5IT,  AIRD,  AISMASS, 
ESERGT,  KE!tGY2.  FlOWl , 
a0«2,  DIGEST,  N0NE2IE3, 
S0HEXG2,  SOOUT.  TESSIN, 
WATT S,  aaCALC 

zlzqi 

ca 

WAT 

Equilibrium  height  of 
liquid  in  the  tank 

TTSNT,  INGEST,  •-•ATLV 

aa 

ca 

GEN 

Height  of  cancer  of 
flow  through  puncture 

tvent,  AIR  LI,  DIEXCT, 
ENESCY2,  flow.,  flow, 

INCEST,  S0NE3E2,  S0SZNG2, 
SOOUT ,  MATIN,  2L2C ALC 

z? 

ca 

GEN 

Puncture  height 

WENT,  AUDI,  ENESCY, 
DIGEST,  NONENEE,  NUOLT, 
VATIN,  ZLHCALC 

a 

ca 

GEN 

Tank  height 

WENT,  DIGEST ,  SOOUT 

zu 

ca 

GZS 

Height  of  water  out¬ 
side  tank 

WENT,  DIGEST ,  SCOUT, 
ZLHCAEC 

r<rc 

ca 

INC 

Height  of  water  in¬ 
side  tank 

WENT,  DICEST,  SCOUT, 
MATIN 

TABLE  C . 6 .  FORMAT  FOR  INPUT  DATA  TO  TVENT  PROGRAM 


Data  may  ba  input  to  the  TVENT  program  through  a  series  of  13  input 
card  sets.  Each  series  corresponds  to  one  data  case.  Multiple  cases  may 
be  tested  by  making  only  one  run  of  TVENT.  To  run  more  than  one  case,  re¬ 
place  the  INPUT  CARD  THIRTEEN  (end-of-file  card)  with  an  end-of -record  card 
(a  7-8-9  punch  in  column  one).  Then  repeat  the  series  of  13  input  card  sets 
for  each  case  to  be  tested.  For  the  last  case  use  the  INPUT  CARD  SET  THIR¬ 
TEEN  instead  of  the  end-of-record  card. 


INPUT  CARD  SET  ONE  (READ  (12,5000)  CID,  CMWT,  CMI,  TC,  CVOL) 


Input  the  cargo  parameters  using  the  standard  FORTRAN  format  state¬ 
ment:  FORMAT  (A3,  7X,  6E10.0) . 

Variable 

Description 

Units 

1.  CID 

Cargo  identification 

— 

2.  CMWT 

Cargo  molecular  weight 

gram/ gram-mole 

3.  CMI 

Initial  mass  of  cargo 

gram 

4.  TC 

Initial  cargo  temperature 

°C 

5 .  CVOL 

INPUT  CARD 

Constant  used  in  calculating  volatility 
of  cargo 

SET  TWO  (READ  (12,5005)  ZT,  VT,  PTI,  HTD) 

Input 

the  tank  parameters  using  the  input  format: 

FORMAT  (3E10.0.I3). 

Variable 

Description 

Units 

1.  ZT 

Tank  height 

cm 

2.  VT 

Tank  volume 

3 

cm 

3.  PTI 

Initial  tank  pressure 

KPa 

4.  HTD 

Heat  transfer  descriptor: 

1  *  Insulated  walls  (ADIABATIC  process) 

-1  *  Isothermal  walls  (ISOTHERMAL  process) 

INPUT  CARD 

SET  THREE  (READ  12,5020)  PDTSC,  ZP,  AP ,  XLP, 

CDORF) 

Input  the  puncture  parameters  using  the  input  format:  FORMAT  (1110, 
4E10.0) . 
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TABLE  C.6.  FORMAT  FOR  INPIT  DATA  TO  TV ENT  PROGRAM 

(Cone 'd) 


Variable 


Description 


Units 


1. 

PDTSC 

Puncture  descriptor: 

1  »  Round  hole 

2  ■  Horizontal  slot 

3  »  Vertical  slot 

2. 

ZP 

Puncture  height 

cm 

3. 

AP 

Puncture  area 

2 

cm 

A. 

XL? 

Width  of  vertical  slots 
(Note:  Enter  zero  for  round 
holes  or  horizontal  slots) 

cm 

5. 

CDORF 

Discharge  coefficient 

— 

INPUT  CARS  SET  FOUR  (READ  12,5030)  SVD,  DPV) 

Input  the  safetv  valve  parameters  using  the  input  format:  FORMAT 
(1110,  1E10.0). 


Variable 

Description 

Units 

1.  SVD 

Safety  valve  descriptor: 

1  »  Unstuck 
-1  ”  Stuck 

2.  DPV 

Safety  valve  relief  setting 

XPa 

INPUT  CARD 

SET  FIVE  (READ  (12,5010)  TINE.  PATM) 

Input 

the  air  properties  using  the  input  format: 

FORMAT  (8E10.0) 

Variable 

Description 

Units 

1.  TETF 

Air  temperature 

°C 

2 .  PATH 

Atmospheric  pressure 

KPa 

Input  card  sets  6  through  11  contain  the  chemical  property  file  con¬ 
stants  used  in  FUNCTION  CPF  and  FUNCTION  CPLHF.  The  variable  T  denotes  the 
temperature  in  ®C  in  the  following  formulas. 


a 
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TABLE  C.6.  FORMAT  FOR  INPUT  DATA  TO  TV ENT  PROGRAM 

(Concl ' d) 

INPUT  CARD  SET  SIX  (READ  (12,5010)  PI,  ?2,  P3) 

Input  the  vapor  pressure  constants  using  the  input  format:  FORMAT 
(3E10.0). 

INPUT  CARD  SET  SEVEN  (READ  (12,5010)  Dl,  D2) 

Input  the  liquid  density  constants  using  Che  input  format:  FORMAT 
(8E10.0). 

INPUT  CARD  SET  EIGHT  (READ  (12,5010)  Zl,  Z2,  Z3) 

Input  the  vapor  compressibility  constants  using  the  input  format: 
FORMAT  (8E10.0). 

INPUT  CARD  SET  NINE  (READ  (12,5010)  Cl,  C2) 

Input  the  specific  heat  of  cargo  constants  using  the  input  format: 
FORMAT  (8E10.0). 

INPUT  CARD  SET  TEN  (READ  (12,5010)  VI,  V2,  V3) 

Input  the  specific  heat  of  vapor  constants  using  the  input  format: 
FORMAT  (3E10.0). 

INPUT  CARD  SET  ELEVEN  (READ  (12,5010)  HI,  H2,  H3,  H4,  H 5,  H6,  H7) 

Input  the  latent  heat  of  cargo  constants  using  the  input  format: 
FORMAT  (8E10.0). 

INPUT  CARD  SET  TWELVE  (READ  (12,5010)  ZW,  VW) 

Input  the  water  characteristics  using  the  input  format:  FORMAT  (8E10 
(8E10.0) . 


Variable 

Description 

Units 

1.  ZW 

Height  of  water  outside  of  tank 

cm 

2.  VW 

Specific  volume  of  water 

cm^/gram 

INPUT  CARD  SET  THIRTEEN  (End-of-file  card) 

This  card  ends  the  data  used  by  TVENT.  Input  a  6-7-S-9  punch  in  card 
column  1  for  the  end-of-file. 


2  '  M  ONI:  NEB 

tiosua  selecis  now  ivpe 

I  -  ILOUl 


miimcriNC  puncvure  an*  iNfouciiN  punciure 


</L(?>.LI.!P>  GO  10 


matatn  ivcni  nm  opio  fin  t.6<osF  tu/ug/u.  ti.gs.or  r»ti 


*00  M.NJ*0-S*<II<2»*1P-»P> 

PINI6-I1. /UI*Ail>*«»l*PA1N**?»  IP*  IN*  I.£-**6*(1W-Il  (*>>•. »/¥«>> 
KI»P=f 


mo, Ian  ivtm  16 in  so/os/n.  n.»>.o/  ran 


dm' tt>  i* »wi'*mi'i*i '»?>»»'«?) am" «?hm'u>  imii  (iio?/'bhhi  i«b  sns 

<  I  »* IS ?B"/ '  HD1VHB01  0// 

in?/'?i )  ii  in 


UCl'D^l  Cl> 

si>  nm 

PA<|-l>*f*A<  n 
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C(**«ijLi»>mxoii)  iNfinodsns 


SUUftOtKINt  CHOKlSf  ?6// 4  OP1-1  MM  H.  13.4i.07  CH( 


NMOI  41 


SUHIOUdtU  CORRECT 


157 


159 


SUAftOUIlME  IUCR6V< I (INI) 


009  01  09(41*  11*0)  UHI 
III  VM  HOI  INlUHll  I IV  3 
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WA*WAO-HV0*Q.0 
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)3N)«jmAN0)  M04  304! 
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^UHMDOllMfc  tHkttfcl  1 6/14  fl*  4 - 6*  4  5  4k  80/08/15.  11.45.0/  ^Abk 
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APPENDIX  D.  SENSITIVITY  ANALYSIS  OF  MODEL 


The  sensitivity  of  the  model  predictions  to  slight  changes 
in  the  initial  conditions  was  determined  by  the  method  described 
in  [2],  That  is,  one  of  the  inputs  supplied  by  the  user,  say, 
the  cargo  temperature,  was  changed  by  some  percentage  and  the 
change  in  the  output,  say  total  mass  released  or  total  discharge 
time,  was  computed.  A  sensitivity  coefficient  was  then  defined 
as 


S  =  (A$/<j>n)  /  (Ax/xn)  .  (D.l) 

where  <p  is  the  output  and  x  is  the  input.  The  subscript  n  re¬ 
fers  to  the  nominal  value  of  x  or  ♦.  A  value  of  S  =  1  means 
that  the  output  changes  by  the  same  percentage  as  the  input, 

S  »  2  means  that  the  output  percentage  change  is  twice  the  in¬ 
put  change,  and  so  on.  Strictly  speaking,  equation  (D.l)  de¬ 
fines  the  "first  order  sensitivity  coefficients"  because  changes 
in  the  output  as  a  result  of  the  interaction  of  changes  in  two 
or  more  inputs  are  not  considered. 

In  this  analysis,  only  the  sensitivity  coefficients  for 
the  input  selected  by  the  user  will  be  determined.  The  sensi¬ 
tivity  to  the  various  internal  parameters  of  the  model,  say, 
the  Chemical  Properties  File  correlations,  will  not  be  deter¬ 
mined  since  these  are  presumed  to  be  "accurate."  Changes  in 
the  numerical  scheme  will  not  be  investigated  either  because  a 
predictor-corrector  type  of  integration  is  employed  in  which 
the  discharge  mass  increment  is  changed  automatically  within 
the  program  to  ensure  that  the  changes  in  tank  pressure  and 
liquid  level  during  each  time  step  remain  relatively  constant; 
this  kind  of  numerical  integration  avoids  the  sensitivity  prob¬ 
lems  of  numerical  integration  that  are  discussed  in  (2] .  The 
sensitivity  of  the  results  to  changes  in  the  various  empirical 
constants  of  the  model  was  discussed  in  the  main  text;  since 
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the  "best"  values  have  already  been  incorporated  in  the  model, 
the  sensitivity  of  the  output  to  these  parameters  is  not  rele¬ 
vant. 

The  sensitivity  coefficients  were  computed  for  a  tank  con¬ 
taining  LNG,  which  is  a  volatile  cargo  displaying  most  of  the 
thermodynamic  interactions  of  interest;  adiabatic  conditions 
were  assumed.  The  nominal  input  conditions  were; 

3  3 

VT,  Tank  volume  =  2  x  10  cm 

ZT,  Tank  height  =  700  cm 

Ap,  Circular  puncture  area  =  2638  cm3 

Zp ,  Puncture  centerline  elevation  =  175  cm 

CD,  Discharge  coefficient  *  0.63 

Relief  valve  =  stuck 

,  Initial  cargo  mass  =  8  x  10^  grams 
T^ ,  Initial  temperature  =  -159. 2 °C  *  114 °K 
patm'  AtmosPkeric  pressure  =  101.3565  KPa 
Tatm'  Air  temPerature  =  21 °C  =  294. 2°K 

All  the  input  parameters  were  varied  by  +5%  except  for  the  tem¬ 
perature,  which  was  varied  by  +1%  of  its  absolute  value;  the 
large  change  of  LNG  vapor  pressure  with  temperature  at  pres¬ 
sures  near  atmospheric  required  a  smaller  percentage  change  for 
than  for  the  other  inputs.  The  relief  valve  was  always  kept 
"stuck"  to  emphasize  two-phase  flow  effects;  it  ought  to  be 
realized,  however,  that  when  the  relief  valve  is  functioning  it 
controls  the  vapor  space  pressure,  and  the  sensitivity  of  the 
output  to  the  initial  temperature  (i.e.,  the  saturation  pres¬ 
sure)  is  much  diminished.  For  all  the  calculations,  the  punc¬ 
ture  was  open  to  the  atmosphere  (nonsubmerged  puncture) .  Note 
that  the  nominal  value  of  temperature  was  taken  as  114 ®R,  rather 
than  as  -159. 2°C,  so  that  the  sensitivity  coefficients  will  be 
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positive  if  the  output  increases  when  the  temperature  increases, 
and  vice  versa. 

Only  the  sensitivity  coefficients  for  the  liquid  outflow 
duration  of  the  discharge  are  presented  here.  After  the  tank 
liquid  level  falls  below  the  puncture,  LNG  vapor  is  slowly  dis¬ 
charged  until  all  the  remaining  liquid  is  boiled  away.  The 
sensitivity  coefficients  for  the  entire  discharge  are,  conse¬ 
quently,  not  very  revealing.  Sensitivity  coefficients  for  the 
total  mass  discharged  are  all  zero  (e.g.,  all  the  initial  mass 
is  discharged,  regardless  of  the  initial  conditions) ,  except 
for  the  two  coefficients  related  to  the  changes  in  the  initial 
mass,,  which  are  both  one.  All  the  sensitivity  coefficients  for 
the  total  discharge  time  are  nearly  zero  (because  the  boil-off 
time  of  the  liquid  dominates  the  calculations) ,  except  for 
those  coefficients  related  to  changes  in  the  puncture  area, 
discharge  coefficient,  and  elevation  of  the  puncture;  these 
are  all  nearly  equal  to  one. 

The  nominal  values  of  the  output  for  the  liquid  discharge 

7 

part  of  the  venting  are;  liquid  mass  discharged  =  5.69  x  10 
grams;  discharge  time  =  82.4  sec;  and  average  discharge  rate  of 
the  liquid  =  6.9  x  10"’.  The  sensitivity  coefficients  for  these 
outputs  are  given  in  Table  D.l,  for  both  increases  (+)  and  de¬ 
creases  (-)  of  the  indicated  input.  A  discussion  of  these  co¬ 
efficients  follows. 

Tank  volume.  Increasing  the  tank  volume  increases  the 
percentage  of  the  liquid  below  the  puncture,  and  similarly  for 
decreasing  the  tank  volume.  The  sensitivity  coefficients  for 
discharge  mass  and  time  therefore  depend  on  where  the  puncture 
is  with  respect  to  the  liquid  surface;  this  explains  the  non¬ 
zero  values  of  the  coefficients.  The  average  flow  rate  changes 
little. 

Tank  height.  Changes  in  tank  height  have  the  same  kind  of 
effects  as  changes  in  tank  volume,  for  the  same  reason,  except 
that  the  direction  of  the  change  of  the  output  is  reversed. 


207 


Puncture  area  and  discharge  coefficient.  These  inputs 
always  occur  in  the  combination  ApCD;  their  individual  effects 
are  identical.  Chances  in  them  do  not  affect  the  total  mass 
discharged.  The  average  flow  rate  and  discharge  time,  however, 
are  almost  related  one-to-one  (i.e.,  S  *  +1)  to  the  changes  in 
Z  or  Cp.  If  the  liquid  did  not  change  phase  and  the  tank 
pressure  remained  constant,  the  coefficients  would  be  exactly 
+1. 

Initial  mass.  Since  the  quantity  of  liquid  below  the 
puncture  is  constant,  changes  in  the  total  initial  liquid  mass 
have  a  slightly  greater  than  one-to-one  effect  on  both  the 
mass  discharged  and  discharge  time. 

Initial  temperature.  Changes  in  the  initial  temperature 
have  a  negligible  effect  on  the  total  mass  that  is  discharged, 
but  have  a  great  effect  on  the  discharge  time.  A  one  percent 
change  in  temperature  changes  the  initial  tank  pressure  by  nine 
percent  and  the  initial  pressure  differential  at  the  puncture 
by  over  fifty  percent,  so  the  large  sensitivity  coefficients 
are  understandable. 

Atmospheric  pressure.  Changing  the  atmospheric  pressure 
changes  the  pressure  differential  from  the  tank  to  the  outside, 
so  the  discharge  rates  and  times  are  affected.  If  air  inges¬ 
tion  had  occurred,  the  sensitivity  coefficients  would  have  been 
larger.  On  the  other  hand,  if  the  vacuum  relief  valve  had  been 
operable,  the  sensitivity  to  atmospheric  pressure  would  have 
been  very  small. 

Atmospheric  temperature.  Since  air  was  not  ingested,  the 
temperature  of  the  atmosphere  had  no  effect. 

It  can  be  concluded  that  only  the  sensitivity  of  the  re¬ 
sults  to  the  initial  temperature  might  create  a  problem  in  prac¬ 
tice.  Even  then,  the  sensitivity  would  be  large  only  for  vola¬ 
tile  liquids  carried  in  pressurized  tanks  or  in  tanks  whose 
relief  valve  is  stuck,  so  that  super-  or  sub-atmospheric  tank 
pressures  are  possible. 
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